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SUMMARY 


This report describes the derivation of the equations of motion 
for a multi-blade rotor. The present analysis advances on 
current capabilities for calculating rotor responses by intro- 
ducing a high twist capability and coupled f latwise-edgewise 
assumed normal modes instead of uncoupled flatwise and edgewise 
assumed normal modes. The torsion mode is uncoupled as before. 
Features inherited from previous work include the support system 
models, consisting of complete helicopters in free flight, or 
grounded flexible supports, arbitrary rotor-induced inflow, and 
arbitrary vertical gust model. These representations are des- 
cribed in previous published texts. 

This report was prepared under NASA AMES Contract Number 
NAS2-6463, "Study, Design, and Fabrication of a High Performance 
Rotor System". Mr. John P. Rabbott, Jr. was NASA Technical 
Manager. 
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DEFINITIONS OF SYMBOLS 

a acceleration, ft/sec^. Also matrix (3.17) 

ajj^ generalized mass in normal mode coordinates, (3.58). 

Also direction cosine (See 7.28, 7.29) 

ap acceleration of a point P relative to the origin 

of the X 5 axis. Fig. 15 

eio reverse flow parameter ,i.:id with steady airfoil 

data, (9.14), (9.15) 

a speed of sound, ft/sec 

ail*» ^ 12 * • 

^ 21 *' ^ 22 * coefficients in expression for correction to 

modal stiffness to account for change in pitch 
(see Sec. 4.2) . 

(aos(l))j^ i-th components of 3-component acceleration vector 

(aoS'^Mi evaluated at O 5 and referred to Xi and X 5 axes, 
respectively, ft/sec ‘ 

a|^ general form of radially independent coefficient 

arising from Ax, Ay, or A 2 components of 
acceleration (See Sec. 7.1c and 7. Id) 

^ blade cross-sectional area. Also, bending 

stiffness matrix (3.10). Also pitch horn 
attachment point. Fig. 24. 

Ax, Ay, A2 linear acceleration terms in acceleration 
components, ft/sec^ (See 7.75 and 7.77). 

A^, A^, Ac ,A^, A^q, A-^, A_^' , A_^, A-^' , A^, 

A-^g, A— -y/' 

orthogonal transformation matrices defined in 
(5.34) to (5.47), and (9.46). Subscripts 
indicate the rotation transformation argument 
angles. 


Ais cyclic pitch control input, rad. Coefficient 

of cosine term 
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'Ji 




B 


%nj 


. ( 1 ) 


®crit 


generalized stiffness in normal mode coordinates, 
(4.9). Also direction cosine (9.89). 

general form of radially independent coefficient 
arising from By, or B, components of 
acceleration (See Sec. 7.1c and 7. Id) 

torsion stiffness matrix, (3.10) 


B^, By, B2 angular acceleration coefficients multiplying 

yiO terms in acceleration components, rad/sec^. 
See 7.75 and 7.77. 




B, . (1), 

limj ' 


R_ . (2) n , . (2) 

' "imj • 




(3) 


B. 

®xm] 


Modal integrals defined in Table 2. 


blade chord, ft. Also, cosine. Also, structural 
damping coefficient 

critical value of structural damping coefficient 


cmn* 

<=d 

Cl 

Cjk 

C 

Cxi Cj 
d 

ds 

ds„ 


general form of radially independent coefficient 
arising from Cx, Cy or components of 
acceleration (See Sec. 7.1c and 7. Id). 

section drag coefficient (See 9 . 11) 

section lift coefficient (See 9.10) 

coefficient in blade modal damping term, (4.70). 

centrifugal stiffness matrix, (3.10) 

C 2 angular acceleration coefficients multiplying z^q 
terms in acceleration components, rad/sec^ (See 
7.75 and 7.77) 

lag damper constant, Ib-ft-sec/rad 

section aerodynamic drag, Ib/ft, (9.11) 

fiber length extending between adjacent beam 
faces after elastic deformation 

fiber length extending between adjacent beam 
faces before elastic deformation 
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D mass matrix, (3.10) 

e blade hinge offset, ft 

e;^ distance of local elastic centroid from local 

elastic axis in the positive yio direction, ft. 

E modulus of elasticity, Ib/ft^ 

E hinge offset vector, (5.20) 

El flatwise blade bending stiffness, Ib-ft^ (6.21) 

Ely edgewise blade bending stiffness, Ib-ft^ (6.22) 

F three component vector of generalized forces in 

loading equilibrium equation, (4.1) . pT » Pj^, 

F2, F3« Also, force acting on a blade section, lb. 

Fl, F 2 » F 3 components of generalized force in force vector in 
loading equilibrium equation, defined by (6.63) 
to (6.65). 

g acceleration of gravity, ft/sec 2 

g* modal damping ratio, = 2c/ccrit 

GJ blade torsional stiffness, Ib-ft^ 

Ti, T 5 unit vectors in xi and xg direction 

Ib blade mass moment of inertia, Ib-ft-sec^ 

Ti'J 5 unit vectors in the yx and directions 

kA structural radius of gyration measured from the 

origin of the yio-ziQ axis (elastic axis) in the 
positive yxo direction, ft, (6.23) 

kyiQr ^zlO radii of gyration measured in the yxo ^nd 

zxQ directions from the origin of the yxo”*10 
axis (elastic axis), ft, (8.13) 

Jcx» k2 1 kyXo» J^zio 
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k 5 unit vectors in the Zj^ and 25 directions 

Kq, Kp, Kp blade leading, flapping and torsion root spring 
constants, £t-lb/rad 

1 section aerodynamic lift, Ib/ft, (9.10) 

m blade mass per unit length sluq/ft (8.13). Also 

aerodynamic moment per unit length, lb. 

quasi-steady section deunping moment per unit length, 
lb, (9.13) 

idq blade reference mass per unit length, slug/ft 

M local Mach numloer for flow perpendicular to 

blade section. Also moment acting on a blade 
section, ft- lb. Also, number of modes. 

Ml, M 2 flatwise and edgewise structural reaction 
moments, ft- lb 


p loading per unit length of blade. Ib/ft 

P distance from the origin of the X 5 axis, ft. 

Fig. 19. Also denotes point of pushrod attach- 
ment (Sec. 12.1). Also intercept of elastic axis 
in x-z plane. Fig. 11. 

q column vector of modal amplitudes, (4.19). Also 

moment loading per unit length, lb. 

qi modal amplitude of i-th iuode. 

Q torsion reaction moment, ft- lb. 

Also position of elastic axis in x-z plane after 
bending displacement. Fig. 11. 


Qj 


Qj 


( 1 ) 


p 


generalized force in j-th modal equation 


( 2 ) 0 .( 3 ) 

contributions to generalized force Qj defined 
by (4.26) to (4.28). Also modal integrals from 
these contributions (Sec. 8 and Appendix 14.3) 


Qj^ aerodynamic contribution to Qj 

Qj^ inertial contribution to Qj 
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Q..A3 


point load (lag damper) contribution to Qj 

aerodynamic contributions to generalized forces 
0 -i'l), 04'2), 0 ^( 3 )^ respectively (See 4.32 to 


4735) . 


inertial contributions to generalized forces 

» respectively (See 4.32 


to 


PI 


4735) . 

Q..P2, Q.P3 


rA 

rp 


pSint load (lag damper) contributions to 

generalized forces 

respectively (See 4.32 to 4.35). 

radial distance in the xs direction from 
hinge, ft (Also denoted by x) 

X 5 coordinate of pitch horn attachment, ft. 

Fig. 24. 

hinge offset, ft. Fig. 1 

X 5 coordinate of point of attaclunent of pushrod, 
ft. Fig. 24 


r^ value of r at blade tip, = R-e, ft 

^ocw value of r at inboard end of counterweight, ft 

R rotor radius * e + rt or * rp + R^, ft. 

R three- component radius vector, (5.20) 

Rl value of xq at blade tip, * R-rp ft. 

Rx, R 2 r R 3 ... modal integrals. Appendix 14.3 
s sine 

Sjk 


coefficient of modal acceleration in modal 
acceleration equation, (4.14) 


Sjk 


( 1 ) 


® jk ^ » Sjk 

contributions to Sjjc from Qj^^, and Qj^^, 

respectively. 
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s 

t 

t 

T 

T 


u* 

U 


Up 

ut 

u 


V 


Vi 


VI 


square matrix of s^t. elements, dimensions 
of M X M, (4.17) 

time , sec . 

forcing function on right-side of modal accelera- 
tion equation, (4.14). 

column vector of t^ elements, dimensions M X 1, 
(4.18). 

steady tension acting on a section of blade, 
lb, (3.5) 


longitudinal displacement ol the blade elastic 
center in the X 5 direction, rt. Fig. 8 . 

magnitude of airflow velocity perpendicular to 
blade section, ft/sec, (9.21). Also when 
subscripted denotes component of total airflow 
velocity, U, ft/sec 

component of U parallel to zio'/ ft/sec. Pig. 18 

component of U parallel to yio'r ft/sec. Fig. 18 

total airflow velocity at a blade section, 
ft/sec, (9.47) and Fig. 19. 

edgewise elastic displacement of elastic axis 
in y direction. Fig. 1, or ys direction. Fig. 8 , 
ft. Also, velocity of a point on the blade, 
ft/ sec. 

v-displacwent of elastic axis in i-th normal 
mode. 

sum of blade elastic edgewise and blade rigid 
body leadirj displacements of elastic center, 
ft, (3.6). 


vii VI displacement of elastic axis in i-th normal 

mode 

(VQs^^Mi velocity vector evaluated at O 5 and referred to 
Xi sucis, ft/sec 

w flatwise elastic displaceneut of elastic eucis 

in z direction. Fig. 1, or :,5 direction. Fig. 8 , ft. 
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Wi 


Wli 

X 


Xo 


Xio 

Xi, X2, 


X 


y 


yp 


yiOcg 

z 

0 ^ 




w-dlsplacement of elastic axis in i-th normal mode. 

sum of blade elastic flatwise and blade rigid body 
flapping displacements of elastic center, ft, (3.6). 

displacement of elastic axis in i-th normal modes 

coordinate. Figs. 1 and 8 . Also vector of dis- 
placements 6 q , wi eind (See 3.10) 

vector of displacements x in i-th mode 

x-coordinate of a point on the blade in the x- 
direction of Fig. 1, neglecting longitudinal 
elastic extension 

vector of coordinates of hub referred to Xj axis 

X 3 resultant torsion moment per unit length, 

resultant flatwise loading per unit length, ard 
resultant edgewise loading per unit length, 
respectively, acting on a blade element. 

position vector with components x, y, amd z. 

Also used to designate force, moment, loading 
an^^ velocity vectors. 

coordinate. Figs. 1 and 8 . Also, vector of root 
stiffnesses, (3.19) 

y 5 coordinate of point of attachment of pushrod, 
ft. Fig. (24) . 

yj^Q value of centroid of section mass, ft. 
coordinate. Figs. 1 and 8 . 

gust inclination, rad. Fig. 20. Also pitch-lag 
coupling angle, rad. ( 12 . 11 ) 

section angle of attack, rad. Fig. 18 

pitch lag coupling angle, rad. 

contribution to isolating terms independent 

of qjj. 
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, o<„(2) , (2) , 0^J3) , (3) 


*rai 


‘m 


mi 


‘m 


‘mi 


Acceleration coefficients defined in Table 2. 


fi 

4 

% 

‘^mn* 



X 


s 


^i 

5ik 

i3- 


blade flapping angle, rad, Fig. 7. Also local pitch 
angle ,rad (in equations for normal mr les of 
vibration) , including built-in twist. 

blade flapping angle displacement in 1-th mode, 
rad. 

lead angle, ^ad. Fig. 1. 

flapping, rad. Fig. 1. 

root end pitch angle, rad. Fig. 1. 

contribution to b^n^ isolating terms independent 
of qk. 

Acceleration coefficients defined in Table 2. 

contribution to c^n^ isolating terms independent 
of qj^. 

blade lead angle, rad. Fig. 7. Also prefix 
denoting perturbation variable. Also elastic 
center d -placement due to torsional displacement 
of a beam element. Fig. 12. 

lead angle displacement in i-th mode, rad. 

Kronecker delta. 

pitch flap coupling angle, rad, (12.11) . 
pitch flap coupling angle, rad, (12.11). 


jv ,Sw elastic center displacements due to torsional 
displacement of a beam element, (5.73), (5.74). 

(S/St)^ differential operator, (7.11). 

^ vector of elastic center displacements in 

twisted coordinates, , A 2 - Also small 

length of blade over which 1^ damper moment is 
applied. 
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^eAcw distance between chordwise location of counter- 

weight and chordwise location of elastic centroid, 
e;^, positive in the direction of yj^Q, ft. 

At time increment, sec. 

Av, Av elastic center displacements accompanying blade 

torsional displacement, ft, (5.75), (5.76). 

Ay* • Az* rotation arms yielding elastic center displace- 
ments due to torsional displacement, (5.67), (5.68), 
Fig. 11. 

contributions to elastic center displacements Av, 
and Aw from the i-th flatwise- edgewise and k-th 
torsion modes, (7.86), (7.88). 


AQ 
Ai, A2 

e 

Cq> ^ 2 ' ^3 
/ 


change in pitch angle 0^ due to flapping and 
flatwise displacement, rad, (12.2). 


elastic center displacements in twisted 
coordinates, ft. Fig. 3 

strain. Also denotes a small quantity. Also 
denotes angle in Fig. 18. 


7 

% 

0 

e' 

~0 

Oo 

^i 


strain coefficients, (6.10). 

coordinate. Fig. 1. Also see Fig. 5. 

coordinate. Fig. 1. Also see Fig. 5. 

y value of centroid of mass of blade section. 

local blade pitch angle due to control input 
and built-in twist, rad, = + 0f ^^9* 

angular coordinate. Fig. 5. Also derivative of $ 
with respect to r. 

total local blade pitch angle, rad, = 

local built-in twist angle, rad. 

sum of built-in twist angle and collective 
pitch control input, '•ad, = 

blade torsion displacement in i-th normal mode. 
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0e 

9a 

0t 

9 


3/4r 




Ai, >2 

/^jk 

/^jk* 

Vmn^ 


I 


• cw 
^ X 

' mn 


(j X 
■'mn 


0 


torsional displacement, rad. Fig. 8 . 
angular coordinate. Fig. 4. 

time-dependent component of pitch angle due to 
cyclic control input, pitch- flap, and pitch- lag 
couplings, ( 12.1 ). 

sum of time- dependent pitch angle 0^ and torsional 
displacement 0^, rad ( = ^ • 

collective pitch input angle measured at r » 3/4R, 
rad. 

pitch angle used to calculate blade normal modes 
of vibration, rad. 

elastic angular displacements, rad. Fig. 8 

correction to modal stiffness accounting for 
pitch difference, (4.46), (4.66). 

contribution to^^jj^, (4.67) 

contribution to isolating terms multiplying 

^k* 

dummy veuriable. Also coordinate Fig. 1. Also 
see Fig. 5. 

position vector with components £, y , J* . 

density, slug/ft^ 

density of counterweight, slug/ft. 

displacement vector 0 ^ 05 , Fig. 15. 

contribution to a^^ isolating terms multiplying 
^k* 

contribution to b|J^J^ isolating terms multiplying 
^k* 

blade torsional displacement, rad. Fig. 1. Also 
inflow angle, rad, (9.20) 

angulau: displacement, rad. Fig. 4. 

angular displacement, rad. Fig. 5. 
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0 ^ 

y'R 

y's 

y' 

00 


blade azimuth angle, rad. Pig. 7. 

angle )^of a reference blade, rad. 

angle of the reference blade at time t » t^, rad. 

angular coordinate, rad. Fig. 4 

angular coordinate, rad. Fig. 5. 

circular frequency of vibration of a normal 
mode, rad/sec. Also emgular velocity, rad/sec. 

components of angular velocity defined by (7.58). 

rotor rotational speed, rad/sec. 


SUBSCRIPTS 


c/4, 3c/4 

j 

k, 1 


e 


m 


n 


quexter chord, three-quarter chord point. 

dxammy subscripts. Denotes direction cosine with 
i and j ranging from 1 to 3. Also denotes mode 
number, with i and j ranging from 1 to M. 

dummy subscripts denoting mode number, with k and 
1 ranging from 1 to M. 

denotes displacement of elastic axis referred 
to ye or ze directions. Fig. 8. Also denotes 
quantity evaluated at hinge. 

subscript ranaing over numbers of elements 
composing Qj®* , , and , Table 2. Also 

first subscript in amn^r bmr> ^> etc. assuming values 
0, 1, or 2 (See Sec. 7.1c and 7. Id). 

second subscript in anm^, b^^^, etc. assuming 
values 0 to 6 (See Sec 7.1c and 7. Id). 


X, y, z, x^, yi, Zi, X 5 , ys, Z 5 , xiq, yio# *10 

Denotes vector component referred to x, y, z, x^,, 
yi, zi, X 5 , ys, Z 5 , xiq, yio» *10 directions 

AIR denotes airflow induced by source independent of 

rotor motion. 

cw counterweight 

H quantity evaluated at rotor hub. 
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I, denotes axis systems. 

Pl> ?2 positions of pushrod attachment, Fig. 24. 

Ox, O 5 quantities evaluated at points Ox and Oc in 

Fig. 15. 


T denotes quantity evaluated at blade tip. 

» >3 , S r € t 6 ^ 0 , , 0* t 0 t ^ 1 1 ^2 ’ ^S ' 

Argument angles in circular functions c and s, or 
transformation matrices. A, rad. 

o denotes a position before elastic deformation. 

Also denotes quantity evaluated at hinge. 

1-10 denotes axis systems. Also 1 and 2 indicate that 

the arguments of matrices A, B, C, and 0 are O 
and respectively. 

SUPERSCRIPTS 


A 


aerodynamic term 


D 


inertial term 


P 


point load (lag deunper) term 


T 


transpose 


X 


superscript in amn^» bmn^» denoting x, y, or 

z (See Sec . 7.1c or 7 . Id) . Also d/d>^. 


D* , Di * , D2* denotes contribution with modal acceleration 
D3* , ( )* terms removed. 


( 0 ! / ( 1 ) , ( 2 ) denotes zeroeth, first, and second order quantity. 
( ) d/dt 


( ) ' d/dr 

(~) complex conjugate. Also denotes vector. Also 

non-dimensional variable. 

equilibrium value. Also running variable 
(Sec. 5.3). 
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abbreviation; J 


e.a. 


elastic axis 


elastic center 
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1. INTRODUCTION 


The work described here is a further step in the evolution 
of methods for calculating the motions and responses of 
rotor systems. This effort advances beyond the single blade 
analysis of Ref. (1) and the multi-blade rotor analyses of 
Ref. (2) and (3), by introducing coupled f latwise-edgowise 
assumed normal modes and a high twist capability in the 
multi-blade analysis. Table A compares the features of 
several rotor response analyses. The present theory 
is incorporated in Sikorsky computer progrcun Y210 (Ref. 5). 
The motivation for the introduction of coupled flatwise- 
edgewise modes is the expectation that a coupled modes 
analysis will be more accurate than an uncoupled modes 
analysis when both use the same number of modes, particularly 
for highly twisted blades. The incorporation of these 
coupled modes with a high twist capability is expected to 
provide a superior tool to handle modo’-n high twist rotors 
and rotors converting to a dual propellor mode of operation. 

Apart from these features, the restrictions and scope of the 
model are similar to those of the previous analyses (Ref. 

1 to 3) . The present model adheres closely to the assump- 
tions of Ref. (1) regarding approximations to terms 
involving products of elastic variables and products of small 
quantities occurring in the equations of motion. These 
assumptions are listed in the next section and are repeated 
in the text as required during the derivations. The proce- 
dure of coupling the rotor with the supporting system was 
inherited from Ref. (2) . The generalized coordinates and 
equations of motions for support systems, consisting of a 
complete helicopter with a rigid fuselage, and a grounded 
flexible support, are described in Ref. 2 and Ref. 3, 
respectively. The method of calculating the normal mode 
shapes and frequencies for coupled flatwise-edgewise -torsion 
motions of a rotating twisted blade is described in Ref. 4. 
This is a finite element (transfer-matrix) method based on 
the Holzer-Myklestad type of treatment. The equations 
were specialized for the present application to yield 
coupled flatwise-edgewise modes and an uncoupled torsion 
mode, as described in Chapter 3. These equations were 
incorporated in the coupled modes module of program Y210. 

Also incorporated in program Y210 is a variable inflow 
calculation module, based on the method of Ref. (5), and the 
gust representation, described in Ref (2) . For completeness 
we have included in this report the equations introducing the 
effects of a cylindrical vertical gust superior to that of 
Ref. (2), but which was not incorporated in the program. 
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F-574 

V-200 

y-141 (Y.201) 

Y-210 

TITLE 

Single Blade Aero- 
elastic Analysis 

Multi -Blade Aero- 
elastic Analysis 

Multi -Blade Aero- 
elastic Analysis 
witn Grounded 
Flexible Support 

Multi -Blade Aero- 
elastic Analysis with 
High Iwist Capability 

B1ad« 

Dynamics 

Externally calculated 
Uncounled Momal Modes 
5 flatwise, 3 edgewise 
2 torsion 

Internally calculated 
Uncoupled Normal Modes 
6 flatwise, 4 edgewise 
2 torsion 

Same as Y-200 

Internally calculated 
Coupled Normal Modes 
7 Modes 

Blade 

Aerodynamics 

2D steady-state airfoil 
data with varying airfoil 
distribution A, B or 

time delay unsteady aero, 
no skewed flow 

20 steady-state airfoil 
data, X, A. B unsteady 
aero, no skewed flow, 
gust input capability 

Same as Y-200 

Same as Y-200 

Airmass 

Dynamics 

Constant or variable 
inflow (calculated 
externally) 

— 

Same as F-574 

Same as F-574 

Constant or variable 
inflow (calculated 
internal .y) 

Airframe 

Dynamics 

None 

Rigid Body Helicopter 
or Compound, Free 
Flight 

Same as Y-200 plus 
ontion for grounded 
flexible «upport 
(10 modes) 

Same as Y-141 (Y-201) 

Airframe 

Aerodynamics 

None 

Steady-State 

Aerodynamics 

Same as Y-200, no 
aero on flexible 
support 

Same as Y-141 (Y-201) 

Control System 
Representation 

Root spring, no 
coupl ing 

Same as F-S74 

Same as F-574 

Same as F-574 

Rotor Control 
Input 

Collective, IP and 2P 
cyclic prescribed huh 
motion or control inputs 

Collective and IP 
cyclic 

Same as Y-200 

Same as Y-200 (Y-231) 

Trim 

Procedure 

Trims to specified rotor 
forces and moments 

Trims to specified 
rotor or aircraft 
forces and moments 

Same as Y-200 

Same as Y-200 

Computer 

UNI VAC 1110 

UNIVAC 1110 
IBM 360 

UNIVAC 1110 
IBM 370 

IBM 370 

Period of Initial 
Development 

1966-1967 

1971-1972 

1972 

1973-1974 

Program Running 
Time {Major 
Iteration S 
Transient 

■ cpu min. on UNI VAC 
1110 

12 CPU min. on 
UNIVAC 1110 

12 cpu mm. on 
IRM 3/0/158 

15 cpu min. on IBH 
170/153 

Rotor 

Configurations 

Articulated, semi- 
articulated, hingeless 

Articulated, semi- 
articulated, 
hingeless, teetering 

Same as V-200 

Same as V-200 

Primary Use 

Blade loads, and 
deflections, pushrod 
loads, blade stability 

Blade loads and 
deflections, pushrod 
loads, rotor hub 
loads, blade stability, 
rotor and aircraft 
gust response 

Same as Y-200 
plus around 
resonance 
stability 

Same as Y-141 plus 
loads and stability 
of high twist blades 




TABLE A - COMPARISON OF FEATURES 
OF ROTOR RESPONSE ANALYSES 
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The present analysis employs Newton's second equation 
to derive the equations of motion. We were inf^luenced by 
Ref. (1) in adopting this method, and by a desire to 
preserve similarities with that derivation, tc assist in 
the development and checking of the equations. The 
resulting equations are transformed into an approximate un- 
coupled set of equations using normal mode coordinates by 
considering the virtual work done by a rotor blade in a 
virtual displacement of a ixsrmal mode. The resulting 
normal mode differential equations are integrated with 
respect to time to obtain the response of a rotor blade. 
Forces and moments induced by the multi-blade rotor are fed 
to the support and the resulting support motion is fed 
back to the rotor to obtain the response of the coupled 
rotor/ support system. 

This report describes the mathematical model in sufficient 
detail to enable a reader to understand, check, and modify 
the theory. The addition of the high twist feature 
together with hub motion terms caused a large increase in 
the number of terms and modal integrals in the equations, 
in comparison with Ref. (1). The proliferation of terms 
emphasized the need for systematic development of the 
equations and exposition of the theory, and the text 
reflects our attempt to satisfy this demand. 

The report is organized on the following lines. Chapters 
3 and 4 give an overview of the approach. It is necessary 
to read these chapters to understand the subsequent text. 

If only general features of the model are of interest, a 
reading of the text may be confined to these two chapters. 

We derive in Chapter 3 the differential equations satisfied 
by the normal modes of vibration and establish the 
orthogonality properties of the modes. These equations and 
their orthogonality properties provide the key to the 
derivation of equations of motion in terms of normal mode 
coordinates from the differental equations for the multi- 
blade rotor. Chapter 4 expands the overview, and introduces 
conventions to describe general terms in normal mode 
differential equations, which determine the organization 
of subsequent chapters, directed to the evaluation of these 
terms in the general equation. 

Chapter 5 introduces the rectangular axis transformations 
employed for the description of the displacement of a point 
on a blade. These transformations consisting of linear 
displacements and rotations, are the generalized coordinates 
in the equations of motion in physical space (in contrast 
to modal space, which is reached after the transformation 
to normal coordinates) . 
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Chapter 6 yields moment and loading equilibrium equations 
similar to those presented in Ref. (1) , to define general 
expressions for the generalized forces in physical space. 
These general expressions guide further the organization 
of the subsequent text, which is dedicated to the systematic 
evaluation of blade element loads in these expressions. 

Chapter 7 begins the evaluation of the terms in the express- 
ions of Chapter 6 for generalized force by obtaining the 
acceleration of a point on the blade, scheduled for sub- 
stitution in the inertia contribution to generalized force. 

In this chapter four different versions of tl.e acceJeration 
are elaborated to ease and systematize the subsequent 
derivation of inertia loads. Chapter 8 completes the 
derivation of the working form for the inertia contribution 
to the generalized force in modal space. 

Chapter 9 yields all the working equations needed to cal- 
culate tlie aerodynamic contribution to the generalized 
force. This chapter includes expressions for the relative 
flow velocity at the blade fur general hub motions and 
gust-induced velocities. 

The text is completed by Chapters 10 to 12 covering the 
following short subjects. These comprise derivations of 

1) lag damper contribution to generalized force; 

2) shears and moments affecting the blade and fixed systems, 
employed to display the response, and required for the 
analysis of the coupled rotor/support system; 

3) expression for contribution to blade pitch induced by 
cyclic control inputs and pitch- flap and pitch- lag 
couplings; 

4) non-dimensionalization rules employed to reduce the 
normal mode differential equations to non-dimensional 
form in the program code. 


3 


PAGE 



\ 


Sikorsky Rircraft 


u 

<^MCRAF 

fi. 


REPORT NO SER-50912 


2. ASSUMPTIONS 


The following major assumptions are made to derive the 
equations of motion. 

1. In general, the rotor hub may be in accelerated motion. 
This state includes unifoxm motion as a special case. 

2. The blade elastic cucis is straight before alastic 
deformation has occurred, and coincides with the feathering 
axis before elascic deformation. 

3. The blade is assumed to have structural and inertial 
symmetry about a major principal axis located on the aero- 
dynamic chord. 

4. Blade flap and lag hinges and feathering bearing are 
coincident. 

5. The blade may have a built-in twist of arbitrary 
magnitude and distribution along the blade span. 

6. The blade angle of impressed collective pitch may have 
an arbitrary magnitude. 

7. Control system flexibility and blade root restraints 
are represented by root springs. Control circuit mass 
and damping effects are net modeled. 

8. Structural damping ratio is the saune in all modes 
of vibration. 

9. Blade element theory describes the aerodynamic forces, 
and radial flow effects are neglected. 

10. In addition to these assumptions, the following 
assumptions are employed to calculate the normal modes 
of free vibration of the blade. 

a) There is no mass unbalance 

b) The elastic centroid coincides with the elastic aocis 

c) Coriolis terms are neglected 

d) Free vibrations occur about a position of ze.ro steady 
displacement. 

11. The following quantities and their derivatives, where 
applicable, eure assumed to be small in comparison to unity. 
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a) Flap and lead angles, in radians 

b) Ratios of translational elastic displacements to rotor 
radius, and torsional displacement in radians 

c) Ratios to rotor radius of chord, center of gravity, 
hinge offser, mass and structural radii of gyration, 
emd elastic centroid 

d) Ratios of hub linear accelerations to /l^R; hub angular 
accelerations so 

e) Ratio of gravity acceleration to^^R 

/) Time- dependent component of pitch angle, in radians 
induced by cyclic control inputs and pitch- flap and 
pitch- lag iouplings. 

12. We list below products of terms neglected as higher 
order terms in the bending equilibrium emd section velocity 
equations. Small quantities considered in these products 
comprise translational elastic displacements, torsional 
displacement, chcrdwlse distances, flap angle, lead angle, 
time-dependent component of pitch angle, hinge offset, 
hub accelerations and angular velocities, and gravity 
acceleration. We neglect in 

a) Flatwise and edgewise bending equations 

(1) Second and higher order products of elastic 
displacanents 

(2) Third and higher order products of small quantities 

b) Torsional equation 

(1) Third and higher order products of elastic 
displacements 

(2) Fourth and higher order products of small quantities 

c) Section velocity equations 

(1) Second 2 uid higher order products of elastic 
displacements 

(2) Third and higher order products of small quantities. 
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3. Differential Equations of Motion for System 
Employed to Calculate Modes of Free vibration 

In this chapter we derive the differential equations 
of motion satisfied by the system employed to calculate 
the modes of free vibration, and we establish properties 
of these equations. These differential equations <md 
their properties nlay a central part in the derivation 
of modal dif ferer .ial equations from the multi-blade 
rotor differential equations. It is recalled that the 
equations used to calculate th^ modes are finite ele- 
ment equations involving tr£uisfer matrices and axe not 
presented as differential equations in Reference (4) . 
Consequently, it is important to verify the forms of 
the differential equations corresponding to the coupled 
modes, emd to establish properties siibsequently used 
to derive modal equations from the differential equa- 
tions for the multi- blade rotor. 

Purposes of this chapter are to: 

1) verify the forms of the differential equation 
satisfied by the system yielding the blade modes of 
vibration , 

2) doncnstrate briefly the property of orthogonality 
which is used subsequently to derive single degree of 
freedom uncoupled mc^al differential equations from 
coupled differential equations, 

3) set out the procedure for deriving uncoupled 
m^al differential equations from the system of coupled 
equations, which is the basis for the later derivation 
of modal equations from the equations of motion for 
the multi-blade rotor system, 

4) indicate the assumptions under which the differential 
equations for the syst«n employed to calculate the modes 
are valid, to indicate restrictions like those required 
to uncouple the torsion motion from the flatwise/ 
edgewise motions, and to show consistency with the 
assumptions underlying the multi-blade rotor differential 
equations , 

5) compare the modal equations of the present method 
with the equations employing the uncoupled modes of 
Reference (1) . We show that our modal equations 
include modes with a predominant ri i body (flap or 
lag) chcuracter, and, hence, that the flap and lag angle 
equations of Reference (1) need not be separately 
derived . 
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Figure 1 illustrates axes systems used for the calcula- 
tion of the coupled modes. The numbering system is 
the same as that used in Reference (4) except that 
primes are added to some subscripted eixes to prevent 
confusion with axis syst^s subsequently used for the 
response analysis. The axes apply to a configuration 
with no steady displacement and coincident hinges. 

Consideration of the equilibrium of forces and moments 
acting on an element of beam yields the following beam 
element equilibrium equations. 


‘ V ' Cc /3 H Pi y + vv ' C- j + cp n a ) 


A 


- O 


{3.1) 






II 


- o 


( 3 . 2 ) 



- (v'f y 4- 4 (% v)'- 



- c 


(3.3) 
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Figure 1. Concluded. 
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structural reaction moments are Q, M 2 illustrated 
in Figure l,e. Applied inertia forces per unit beam 
length aire p^r Py# Pz. Applied inertia moments per 
f unit beam length'^are qx# qy> qz* Steady quantities 

are denoted with a hat, ^ma perturbations are prefixed 
with S[Spx ^Pyr etc.), or are untreated (u, v, w) . The lo- 
cal pitch angle is /£ and is not to ^>e ronfused with the 
flap angle of the multi-blade rotor response equations. 

1 ^ These equations could have been derived also from tJr,e 

beam element equilibrium equations of Houbolt and Brooks, 
equations (17) and (11) , Reference (6) , by dividing 
quanticie& into steady and per tuba tion values. 

The structvtral reaction moments may be expressed in 
terms of section stiffness properties and displacement 
derivatives by means of equations (15) , (16) , and (17) 
of Reference (6) . 

The inertia loads Px» Py» Pz» qx» qy» qz derive from 





(3.4) 


In these expressions, y and z are displacements of 
an elemental mass in a blade section after elastic 
displacements occur, and ax, ay, az are acceleration 
components of the elemental mass. Neglecting blade 
longitudinal extension u, the accelerations apply to 
a point distance Xq from the origin of the x-y-z axes. 
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Figure 2 illustrates inertia force and moment positive 
conventions, and location of the point where forces, 

Px* Py' Pz and moments, q„, qy, are evaluated. This 
point' is on the x-eucis and is^not at the elastic center 
of the blade which, in general, is displaced to some 
other point. The purpose of the illustration is to 
emphasize that the applied forces and moments of the 
coupled modes equations are defined differently from 
the forces and moments of the response analysis, also 
designated px, Py, Pz, qx> qy> qz» which are evaluated at 
the elastic center and resolved to the 5 system of 
the response equations. Also, qy in the response 
equations is opposite in sense to the qy of the coupled 
modes equation. 

Employing the assumptions listed in Chapter 2, we 
derive below from Reference (4) results for inertia 
forces emd moments. Included in these assumptions cure 
the following assumptions consistent with those of the 
response equations. 

1) Fourth order products of elastic displacements and 
small quantities, like radii of gyrations, are neglected 
in the torsion equation 

2) Third order products of elastic displacements and 
small quantities are neglected in the flatwise and 
edgewise equations. 

In addition, the following assumptions are made to 
uncouple the torsion equation from the flatwise - 
edgewise equations, and to produce real elgensolutions: 

3) - 0, no mass unbalance 

4) Sa " 0, elastic centroid coincident with elastic 

axis 

5) Coriolis terms are not included in the inertia 

forces. This is an assumption of 
Reference (4). 

6) Free vibrations occur about a position of zero 

steady displacement. 
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The following results are obtained: 


-A 
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- o 


P\i O 
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SQ -- (& J + T K^) 
r + r-p ) 

>5 -o 
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lx -■ 




O 


r O 


(3.5) 


% IP z O 

-'"’‘oPc. 

C ; - »n A'”(xo+ Tp ) - xnw 

fz 

Sq - - >^0 + V 
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The coupled nodes nodule can embrace less restrictive 
models than those satisfying assumptions 1) to 6) . To 
prevent the use of an inconsistent model, inputs like 
^t' ^z' ^^1 defined in Reference 

( 4 ) , are overridden by instructions in the program 
setting these parameters to zero. This has the effect 
of yielding (3.5) exactly. The differential equations 
satisfied by the coupled modes are then exactly those 
derived in this report (see 3.7 to 3.9, below). This 
observation establishes the consistency of the equa- 
tions for the natural modes in the progreun and this 
report. 

Substitution of (3.5) in (3.1) to (3.3) yields the 
differential equations of motion satisfied by the 
system employed for the coupled modes analysis, ouid 
the substitution 

(3.6) 

P’e 

W/ ' w p j= ^ o 

V, ^ 


transforms the equations into 

(3.7) 


- Sv, = 


-((G3^T^p=) e/)' 
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1 1 


(w/t) 


•f K>i W 


- o 


(3.8) 


•SX, . 


-(v/f)'- 




/ ( 


O 


^Jl i hr\ 


(3.9) 

These equations aire expressed in iaatri> form, to effect 
the most convenient derivation subsequently of the proper- 
ties of the system. 


(3.10) 


(A.")"+(B.')'+Cx f Oy .-O 

A s O C> C> 

^ ^ ^ ^ / 3 ^ /3 ) ^ /3 

B Z o O 
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-k: o 
O o 

o o 

O 

o ^ 

e;,w, , 


It is seen from (3.7) to (3.9), or (3.10), and 
bound 2 u:y conditions (3.15) and (3.16), Section 3.1, 
that the torsion motion is confined to (3.7) which 
does not involve vertical or inplane motions, indicating 
an uncoupled torsion motion. In general, for non- 
vanishing pitch,/^, vertical and inplane motion terms 
appear simultaneously in (3.8) and (3.9) indicating 
that these are coupled. The uncoupled character of 
the torsion motion euid the coupled character of the 
vertical and inpleme motions are consequences of 
assumptions 1) to 6) , listed above. Equation (3.10) 
does not indicate fully coupled motions and should be 
viewed as a convenience. Orthogonality results 
subsequently derived from (3.10), combining torsion and 
other displacement modes in one expression, are not 
to be construed as being for fully coupled motions, 
and apply to uncoupled torsion modes, as justified 
above. (The orthogonality condition can be bro)cen up 
into terms only involving torsion modes, and terms 
only involving vertical and inpleme modes, to indicate 
explicitly the uncoupled character of tiie torsion mode. 
This is not done in the large body of text because 
no advcmtage derives from doing so for the present 
exposition. Computational efficiency could be improved, 
on the other hand, by distinguishing the uncoupled 
character of the torsion mode, as discu»<sed sub- 
sequently . ) 
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3. 1 Boundary Conditiong 

Equation (3.5) yialds displacement derivatives as 
functions of structural moments. 

- 5c?/rfi3 4 ) 



_j 




Ll 1 ^ Sj3 




(3.11) 

(3.12) 


At the hinge, including root springs, we have 



(3.13) 


rc^ 



^ K, /3f ^ 

(3.14) 

L ^/3 

- 


's 5 



Substitution of (3.13) and (3.14) in (3.11) and (3.12), 
respectively, yields 



^/o J [(3 J ~h T 



(3.15) 



(3.16) 


where a is the matrix 
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(bi 


^ ^ Lrl 6 ^ . ei )s o 




V , £l, 


For subsequent substitution of the boundary conditions, we 
note that A of (3.10) also is 


o > o o 


Eiuploying (3.16) we obtain 




0 \ O O o \ o o 


O ' CL 


^0 ftoJ 




f " 

(^oFd 


This is rewritten as 


( ")c 


(3.19) 
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with 

b - O , Po 


At the tip, s • Sm 2 0 and from (3.11) emd (3.12), 
we find w" ■ v" ■ 0. Also, shears, «* SUy' “ 0 

at the tip. Differentiation of (3.12) then yields w ■> 
V a 0, at the tip. Using (3.6) we can state the 
boundary conditions in terms of 0^, Wj^, Vj^. 


xq a 0 (hinge) 


(3.20) 

6e -- /3p 


--0 

I ) I C Cl ; f^jr ^ P>Q 



K 1 

K^l^o) 


with a defined by (3.17). From this last equation we obtain 



with 

(^"*1 o j kp /3p ) ‘^o 

and A defined by (3.10). 

Xq ° R l. (^ j 





w, 


1 1 1 


z V 


I 1 1 
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3.2 Orthogonality of Modes 

For a system in free vibration at frequency (3.10) 
becomes 




/)/•) " i-C&x'Y + C)f- 


o 


(3.21) 


To establish orthogonality of modes we evaluate the 
work done by a mode with f requency , in a virtual 

displacoaent xj. To evaluate the work qucmtity, we 
premultiply (3.21) by x? and integrate from hinge to 
tip. We obtain 




i '[ f x/ Ce ) V x^f cx, -co,^x/o>c^ J 




(3.22) 


Integration by parts twice of the first term yields 

Invoking boundary conditions (3.20) we are left with 

"\'‘l f r^l -r" » 

K ^ 4 ^ 4 -Xj P. X^- dK^ 


Again, from (3.20) 
and (3.24) becomes 


(3.23) 


(3.24) 


( 3 . 25 ) 
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// 



(3.26) 


Integration ny parts of the second term in (3.22) 
yields 






(3.27) 


j^n view of boundary conditions (3.20) emd the fact that 
T = 0 at Xq = we find B = 0 at Xq * Ri and (3.27) 
becomes 


X )o~^c '^J o 


(3.28) 


The boundary conditions (3.20) yield 


Kpp^^pp. 


J 

(3.29) 


and (3.28) becomes 


•L ^ ~ ' "o 


'^1 -r'- ' 

-X- 0>Tcp (Sif 


o J 


o 

(3.30) 


Equation (3.21) for the virtual work becomes 

<p Pfo- + ^Fi ^ I'^o^ Po 

(*? 




'J 


0 .■ 

J . t ^ I ~ ^ V ®X.- -f- 




(3.31) 
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Sixailarly the work done by the forces of the j-th mode 
in the i-th virtual displac^ent is 


Kp/3^^. /3^-. h Kf/3p_./3f_^ t ^oP‘d -^13 

L 


D 


B xj f P 


y 

(3.32) 


where ^2 is the frequency of the j-th Dode. At this 
point we note that products like x^”AxT" eure scalars. 
The value of the product is unchanged By taking its 
transpose. Thus 


I '/ us. 

A X ' J 

^ j 


T 


T \ T " 


Since matrices A, B, C, D are symmetric 


A ^ s'" c*" 


A , 13, C ^ D 


(3.33) 


(3.34) 


and it follows 


r '' 






A Xj J - , f\ 


-X, 


(3.35) 


Taking the transpose of (3.32) and using (3.35) and 
similarly treating other terms in (3.32) and sub- 
tracting the transpose of (3.32) from (3.31) we are 
left with 


( 




UJ 




0 ?c / ' O 


(3.36) 


If eigenvalues and eigenvectors are real and distinct 
this requires 






o 


^ +J 
^ ^ J 


(3.37) 


22 

PAGE 


$A 29 REV 0 


I 


1 


Sikorsky 


Pircraft 


u 

0. 


REPORT NO SER-50912 


This is the statement of the orthogonality condition 
of the (normal) modes and it applies to the coupled 
system (torsion is uncoupled) with or without root 
springs. Expanding (3.37) we find 



f- w/, W, . ^ Vj , 
«. J 



Z o 


(3.38) 


For the rotated system of axes of Figure 3 defined by 
the transformation 




a -- 

1 

C' 


Ae ) 

U/, 

o 

Cp 

S/3 

K \ A, f 

1 

0 

- ^/3 

o 

1 

\ A ^ / 


(3.39) 



A 


(3.40) 


(3.41) 



\ o o 

O C/3 

O " ^/3 


(3.42) 
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Figure 3. Rotated Axes Defining Elastic Center 

Displacement in Coupled Modes Analysis. 
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we find 





^ oA ^ 

o 


(3.43) 


Equation (3.37) becomes 


^ x^- ' - X A = 

r d <- ^ j 

/ 0 t -- j 


(3.44) 


Equation (3.44) shows that the form of orthogonality 
condition (3.38) is unchanged under rotation. Equation 

(3.44) is the form used in the coupled inodes module. 


3.3 Reality of Eigenvalues and Eigenvectors 

The assumption that eigenvectors X and eigenvalues u / are 
real is proved below. It is recalled that this assumption 
was used to prove orthogonality of modes. Denote A > -or. 
Denote complex conjugates by an overbar. Equation (3.21) 
yields the following forms 

) -i (B X ' ) -fCxfAO-x so 

-^( 6 ^')' ^ C} f A Ox -- O 


(3.45) 

(3.46) 
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Premultiply the first equation by x and integrate from 
0 to We obtain 

/ xYOx') ^ ^ Ax Yy Jo^y^-O (3.47) 


Invoking boundary conditions as we did to derive (3.31) 
and noting that boundaury conditions are real, we find 


f^Pj * ^ 1^‘^J 

X X ^ A X 5 y ~{ X Cx -f A > OyJ cS'Xq 


O 


(3.48) 


Similarly a premultiplication of (3.46) by x*^ yields 

Kf, Prc Pfj ^ Pf^ Pfj + ^oPocPo^ 

'ox' -t ° 

O ^ -1 


(3.49) 


In a similar way to the derivation of (3.35) we find 



(3.50) 


Taking the transpose of the terms in (3.49), using (3.50) 
and similarly treating other terms, and subtracting the 
transpose of (3.49) from (3.48), we obtain 


(A -a) 


r"' 


- r 
X 


Ox ci> 


o 


(3.51) 
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The integral is not zero in general and consequently 

^ (3.52) 

which proves the reality of the eigenvalues. Since 
matrices A, B, C, D are real in (3.45), the substitu- 
tion of a real eigenvalue in (3.45) can only yield real 
eigenvectors. This shows that eigenvectors are real. 


3.4 Derivation of the Modal Equation 

In this section we illustrate the procedure for deriving 
the modal equations, which cure single degree of freedom 
uncoupled differential equations, from t^le set of coupled 
differential equations (3.10), by invoicing orthogonality 
cxsndition (3.37). The procedure is the basis of the 
later treatment of the response equations whose purpose 
is to generate approximately uncoupled modal equations 
cund thereby facilitate and stabilize the integration 
with time of system response. 

Express the displacement in (3.10) as a sum over the 
number of modes of modal vectors Xj^ multiplied by modal 
amplitudes qj,. 






(3.53) 


It is cx>nvenient to use the summation convention to 
contract the algebra. 

(3.54) 

The appearance of a repeated subscript will always imply 
summation on that subscript in this report and this 
convention is extensively employed. We recall that here 
xi is a 3 component vector ^md qj^ is a scalar. 

T 

Substitute (3.54) in (3.10, premultiply by Xj and 
integrate from 0 to Ri. This yields 
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Xj D ^ 0 ' ^ 


(3.55) 


Repeat these steps with equation (3.21). This yields 


/ ^7CB-;c7)V 


-t 


T 


Ox.y^cily^r O 


X ifcj dl. 


(3.56) 


Subtracting (3.56) from (3.55) yields 




f. ('^C ) - 


O 


or with 


(3.57) 




~ f X-Dy:^<^^ 

O J <- o 


(3.58) 


we obtain 


CljC Ji + ^ ^“-JC 


z O 


(3.59) 


Since 


^ J -■ 


(3.60) 


in view of (3.37), equation (3.59) is the desired uncoupled 
single degree of freedom modal equation. 
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If equation (3.10) had been 

+ (6^-3'4C^ + 0;^- -- F 


(3.61) 


where F is a 3 component vector representing additional 
forces, like aerodynamic forces, the modal equation wou' i 
become 





f to Q . 





O 



(3.62) 

(3.63) 


Our approach to the equations to be used to find the 
response of the rotor is based on the objective of creating 
a form very similar to (3.62). This is the guiding 
principle of the subsequent derivation. After modal 
responses are known, physical displacements, velocities, 
and accelerations derive from equations similar to (3.54). 


3.5 Cgnparison of Modal Equations 

We show here that the set of modal equations of the subject 
method includes the flap and lead angle equations of 
Reference (1) . Consequently, it is not necessary to 
derive separate equations expressing the equilibrium 
of hinge moments (the flap and lead angle equations of 
Reference (1) ) , as was done in Reference (1) . 

Consider a (first) mode which is essentially a rigid body 
mode. The eigenvector of modal displacements is 


- c;. r- 


o 


0 


(3.64) 

(3.65) 


following normalization to flap eingle is 




(3.66) 
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(repeated suffix liidicates siunmatlon on i over the modes) . 
Flapping participation is small in the other modes emd it 
follows 



(3.67) 


and substitution of (3.65) and (3.67) in (3.62) yields 



(3.68) 

(3.69) 


= blade moment of inertia. 

This equation essentially duplicates the flap angle equa- 
tion, 1-91 of Reference (1). It may be shown similarly 
that the lead angle equation, 1-96 of Reference (1) , 
is essentially duplicated by a member of (3.62). Thus it is 
not necessary to derive separate equations expressing the 
equilibrium of hinge ’loments, corresponding to the flap 
and lead euxgle equations of Reference (1) . 
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4. General Features of Modal Equation for Multi-Blade 
Rotor System 

We discuss here general features of the blade modal 
equation applicable to the multi-blade rotor system, both 
to summarize our approach and to organize subsequent 
chapters, which are directed to the evaluation of terns 
in this general modal equation. 

The differential equation ratlsfled by the motion of each 
blade of the multi-blade rotor is assumed to be 

( X " ) ^ ( 5, X ' ) + C, X t 0, X ; F ,4.1, 

This equation may be verified for the multi-blade rotor 
by isolation of terms on the left- side from the loading 
equilibrium equation established in subsequent chapters. 

The terms on the left-side &xe first order structural. 
Inertial, and centrifugal terms that would occur also 
for a motion of the system in free vibration, considered 
in Chapter 3, and presented in (3.10). As a result, 
these turms will exhibit the desired orthogonality 
properties established for the modes of free vibration, 
enabling us to form a modal equation like (3.62). The 
terms on the right side of (4.1) contain all terms except 
those on the left-side. These terms include aerodynamic 
and concentrated applied loads, like those due to the lag 
damper and inertial loads , including Coriolis terms and 
terms arising from mass unbalance and non-coincidence of 
elastic centroid and elastic axis. 

As a preliminary to the derivation of the equation of 
motion in modal coordinates, we state that a mode of 
free vibration satisfies 

( Aj X '')"*( Oa = o 

which is (3.21). 

The different Rubscripts on the matrices in (4.1) and (4.2) 
distinguish different pitch Emgles. Subscript 1 signifies 
that the argument angle in Ai, , C^ and (B^ and 
are not functions of pitch angle; is 
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© -■ e, , 

with 

<9c " K * ®e 


(4.3) 


(4.4) 


The angle Is the resultant local pitch, summed from 
the collective pitch input, (? 3 / 4 r, and built-in twist, ^ 3 . 
The angle 0^ is a time dependent contribution to pitch 
angle arising from cyclic control input, and pitch- flap 
and pitch- lag couplings, such as will be defined in a 
subsequent chapter. 


The subscript 2 in (4.2) signifies that the argument 
angle in matrices A^, B 2 , C 2 , and D 2 , is the pitch angle 
used to find the modes of free vibration, which we shall 
designate ^ 3 *. (The angle in (3.21) ). 

The distinction between the pitch distribution used to 
find the response of the multi-blade rotor and the pitch 
distribution used to find the modes of fre^ vibration is 
necessary because these distributions are different, in 
general. When cyclic inputs and pitch-flap and pitch- lag 
couplings occur 0 cannot be the same as 0^* owing to the 
fact that 9 embodies a time-dependent component 0^ 
while 9n* is time- invariant. Even in the case of no 
cyclic Inputs and no pitch couplings, a requirement^* 3 9 
would lead to a recall of the natural frequency calcula- 
tion module wnenever collective angle ^ 3 / 4 ^ ch£uiges (such 
as in a "major iteration"), adversely Increasing program 
execution time. To avoid the calculation of normal modes 
whenever ^ 3/43 is changed, and to allow for time-dependent 
differences' between 0 and $^* , the approach will be to 
distinguish between 0 and ^ 3 * and properly correct for 
the effects of these differences. It will be seen that 
the nature of the correction is an addition to the stiff- 
ness (or frequency) of the natural mode, accounting for 
pitch differences. It is assumed that the normal mode 
shapes are the same at 0 and . This assumpt j on is 
embodies in (4.5), given immediately below. 

Proceeding now to derive the modal equation, we follow 
the step.i of Section (3.4). We express the displacement 
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e' 

X-- X.?<: 


(4.S) 


(see equations 3.53) and (3.54) ), premulfply (4.1) 
and (4.2) by x-*^, integrate their terms from 0 to r^ 
and subtract tne resulting integrals. We obtain 


“O' C + / V " \ " 


(4.6) 


after substituting Bx » B 2 » = D 2 * Equation (4.6) is 

re-written 




(4.7) 


( D, \ 

(4.8) 


3 

- L/<^ 1 

J c, 

/• r' -r / 1' \ 

(4.9) 


' •[, 1^-’ ^ 



+ ■Yj 

(4.10) 

(\ 

C't t- _ . 



- j ^ ^ 

•^0 ^ 

(4.11) 
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Equation (4v7) is the same as (3.62) except for the 
occurrence of the additional stiffness terra^ji. In 
view of the orthogonality property of ajj^, described in 
(3.60), we can aver that with^ji and Qj excluded, equation 
(4.7> would be an uncoupled single degree of freedom modal 
equation. The termsy^ji amd Qj couple the set of equations 
(4.7). An important assumption for the success of the 
time integration of the motion is that these couplings 
are always sufficiently weak to permit the generalized 
mass, ajj^, and stiffness, bjj^, to dominate the motion in 
.aodal coordinates, q^, and he'’- ;e to preserve essentially 
the uncoupled character of (4.7). 

The interpretation ot/U^^ is that this is a correction to 
the stiffness bji of the coupled modes of vibration, which 
accounts for differences in structural and centrifugal 
stiffnesses corresponding^ to the actual pitch 0 and that 
used to find the modes, These interpretations are 

amplified in a subsequent section of this chapter, where 
is also reduced to working form. 

To use (4.7) in our integration of the motion, we turn it 
into an equation in qi, q^, and qj^ by expressing in Qj 
displacements, velocities, £uid accelerations as modal*^ 
sums. Typically, flap angle is and vertical 

elastic displacement is w = w^q^. Next, (4.7) is 
re-arranged to place modal acceleration, q'j^, on one side 
and mode' cunplitude, qj^, and velocities, qj^, on the other 
side, and this form combined with time- integrations yields 
the history of motion. 

To effect a form of (4.7) with all terms involving q on 
one side, generalized force Qj is broken into inertial, 
aerodynamic, and point loads ^mponents. The point loads 
are loads induced, for example, by the lag damper. With 
D, A, P denoting inertial, aerodynamic and point load 
forces this yields 

The only additional source of modal acceleration terms, 
explicitly occurring in the modal equation, is Qj^» 

For example, terms involving g can come from mass 
unbalance. The residue of Qj° it is deprived of 

acceleration terms is wriiten Qj^ . We then express 
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0 


o 

<5 , 


( 0 . 
jk'j 


S 


JK 1 


U) 


K 






(4.13) 


where are acceleration coefficients 

deriving from torsion, flatwise, and edgewise equations 
respectively, explained subsequently in this chapter. 
Replacing the dummy subscript i by k in (4.7), we re-express 
(4.7) as 




J 

(4.14) 



s s s 

7 J K ^ K 

(4.15) 




Ik 




(4.16) 

Letting 

M 4= number of modes, and 



,? - 

r S' - ] 

^ ^ Kj -Hx H 

(4.17) 



i 5 -M X 1 

(4.18) 


^-1 ^ 

1 ! -Mx ( 

(4.19) 


the matrix equivalent to (4.14) reduces to 

Z 5 ' V ^ ^ ' (4.20) 


and this is the desired form for the modal acceleration. In 
keeping with the expected essential uncoupled character of 
(4.7), matrix S, although fully populated, will be 
diagonally dominant. 
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Modal velocities and amplitudes derive respectively from 

^ ^ ^ - j^(t) I (4.21) 

"|c(-t+At) : ^ -^7< (4.22) 

User experience with rotor response equations employing 
uncoupled modes suggests that (4.22) is superior to a 
Taylor series, and this is the justification for its use 
here. 

Equation (4.20) is the target of the subsequent text which 
seeks to derive S and T in (4.20), and which is organized 
in an orderly manner to do this. 


4.1 General Expressions for Generalized Forces 

To aid the rational organization of the work required to 
form the elements of the modal acceleration equation, 
(4.20), it is of value to break the generalized force 
into terms of distinct type, namely inertial, aerodynamic, 
and concentrated load. This section lists the conventions 
used subsequently to define the components of generalized 
force. 


The generalized force in the modal space, Qj, is 





-X • F cil'^ 

J 



f V, • 




F dL 


(4.23) 

(4.24) 

(4.25) 

(4.26) 
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(4.27) 


(4.28) 


Element oi is the work done by torsion moments (per unit 
length of ''blade) moving through the torsional component 
of virtual displacement, (9j, of the j-th mode. Similetrly 
0 j 2 ) and are the works done by vertical and inplane 

forces per"* unit length moving through vertical and inplane 
disDlacements, j , and of the j-th mode. Element 

may be thought of as a generalized torsion moment 
applied to the entire blade, and Qj 2 ) and Qj * are generalized 
hinge moments. It is recalled that the forces F]_, P 2 , F 3 
embody inertial forces in addition to aerodyn 2 unic and 
concentrated loads forces, and these forces are not to be 
construed as only of external origin. 


Decompositicn of forces into iiiertial, aerodyneunic and 
concentrated load components is effected by the following 
definitions. 



t ^ 1 

P 

f F 
1 

(4.29) 


/- 0 r A 

^ ^ 

f/ 

(4.30) 


^ D ^ A 


(4.31) 






Superscripts D, A, and P designate terms of inertial, 
aerodynamic, and concentrated load origin, respectively. 
The generalized force components may then be written 



to 


(i) 






J 







(4.32) 

(4.33) 
(4. 34) 
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where 




^ f , 0 
-I J. V,-; / , , 


C ’M 


J ' ^'z 


J 3 


<5.“' 

<5,'’’ -/;v f u 


l; 'j ^ 




o J 


'^/-''/'V F/A 

o 'j ^ 

op-rvF.''* 
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The Interpretation of a typical term is that Q 


,d2 


is the contribution to the generalized force by the 
vertical inertial loads, ^2^* moving through the vertical 
dlscplacement component, w^j, of the j-th mode. The term 
is a similar contribution deriving from aerodynamic 
forces. 


In view of (4.12), (4.25) and (4.32) to (4.34) we may 
express inertial, aerodynamic, and point loads contribu- 
tions to generalized force Qj as sums of contributions 
from torsion, vertical, and inpleme forces. 





J 

(4.36) 


(4.37) 

4 C- * Lj 

(4.38) 


The residues of Qj , Qj , and Qj*^-* after they are 
deprived of modal accelerations are written Qj“^ , Qj^^ » 
and Qj°^*, and the inertial generalized forces may be 


written 


(('} '■'i 

S' 


1^, 

o. 


- u 






Q • 

^ J J 


- 5 

jK 

Ca),. 

- V/t Ik 
5 

■ j'K 


(4.39) 

(4.40) 

(4.41) 


It is seen from these definitions that si]^ , Sj]^ , and 
sj^' are acceleration coefficients whose*^ sources are 
torsion, vertical, and inplane forces, Pi*^, F 2 ®, F 3 ®, 
respectively. The repeated suffix k indicates summation 
on the number of modes. Defining 







(4.42) 
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and substituting (4.39) to (4.41) in (4.36), we obtain 


0 , 


0 


c 


D 


5 .' 


( 0 . 


'J J k 

which is (4.13). 






(4.43) 






Because the forces F 2 » and F 3 contain a large number 
of terms, which will appear in subsequent development, 
reduction of execution time becomes a worthwhile goal. 
Such a reduction can be achieved by distinguishing the 
uncoupled character of the torsion mode, in the formation 
of generalized force Qj. Recognizing that 



(4.44) 


for a torsion mode, and 







(4.45) 


for a coupled f latwlse-edgewise mode, we indicate a means 
for effecting a reduction of computational effort needed 
to form Q j . 


4.2 Correction to Modal Frequency to Account for Pitch 
Differences 


In Section 4.0 we showed that a consistent derivation of 
the modal equation, employing coupled assumed modes 
calculated for a pitch distribution 0^, gave rise to a 
correction to the modal stiffness bj}^ to account for 

the change in"^ pitch to $ during the motion. This correction 
may be viewed also as a modal frequency correction. This 
section describes the derivation of the working form of^</j]^. 

The expression for^i/jv, equation (4.10), is the starting 
point for the derivation of the working form. 
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Using identities Ely = EI^, EI^ = El 2' ^zlO “ ^2' ^ylO “ ^1' 
we obtain from (3.1u) 


(4.47) 


hl.C^ ' t I 2 s, 


O t 1 Q 


(t 




(4.48) 


O 

fc-Xcj + e I2 ^ V 


(£I ■ £'Iz)\. ^ 






(4.49) 


w.n,"c /■ 


^ t o ^ “*e 


W\ J^L C - 




^Sb" 


- V-r^Ti 


(4.50) 
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By means of (3.24) we can show 






Of, 


"4 , 


From (3.19) 

^ k ,/3 


o 


r* '^ 0 ' '' A V 

j (p> 

^ r\ 


ciu. 




^ n S 


O ,^frK ) D"^K 


K 

(4.51) 


(4.52) 


(4.53) 


where >^Emd S signify flap and lag angles (replacing the 
notation j5f and of Chapter 3). Equation (4.51) becomes 



T 



tj* ocj ’'ffl,- 

•/ Q 


(4.54) 


With 






(4.55) 


equation (4.54) reduces to 
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o 




+ tf h * '^-9 <- J 9 ^/ 


with 


( 4 . 56 ) 


(4.57) 


a 


-)( 


I 




^<9^* 


<X, 


I 9^ 


eljflj 


-(ei, 

(lI^ -(f ^ 5 ^ C^(E1 . 5^ ^ ^ tl^ ^ (9 x) 


-(ei 


j ^9 ^ e 


^/3 

^ r= o 

)(ei,-El;^)%^^Cp 


/3 


a I 


‘I 




yz O 

:(ei^c= . £ I ^5 = )( £1, ■ e i“5 ) * 

' (£li-£i'j) 5©'-© C t I ^ Cgjj K *. 1 12 ) / 

1 2^ ■ I ^^,3’' *^ee ^ 


(X ^ 


EX 2 1 i> 


0 ^ ^ 

^6 "X ■(■ 

^6 


I 

- I 
r ^ O 
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The integral in (4.56) is 


f* '*'T -r- \ l< 






y <9 ' y ^''(9 


W. 


i! 


K 




The second term in (4.46) is 


J nr^Jl <},J{ 3s *^©/3 / J ^ 


(4.58) 


(4.59) 


Substitution of (4.56) to (4.59) in (4.46) yields the 
expression for^jj^. 
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From a computational viewpoint, (4.58) and (4.59) are not 
efficient because the presence of a time-dependent part 
in 0 i = Of, + fij;) would require these blade integrations 
to be performed at each increment of time of the motion. 
To permit employment of a more efficient calculation of 
integrals are developed below, independent of time, and 
which need to be evaluated only once before the motion is 
calculated. 

Expressing 


6' c'c + 


(4.60) 


we assume 0^ to be small, specifically of the order of 
the small quantities, li)ce c.g. offset, and section radii 
of gyration, and we neglect high order terms contributed 
by/Ujj^qj^ to the modal equation of motion. To exploit 
the different orders of approximations assumed for the 
torsion loading equilibrium equation in comparison v/ith 
the f latwise-edgewise loading equilibrium equations, we 
recognize that the uncoupled character of the torsion mode 
transforras these approximations into equivalent 
approximations in the modal space. In the torsion modal 
equation 


G' 


(^j.< * 

^ 


0 

K ■- ^ K 


hjK m J )(c 

O / o ' ^ 


- C -5 ^ -X 


ais 
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Cl 


{<) 


V I ! 


O 'J 


(4.61) 

(4.62) 

(4.63) 

(4.64) 
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we neglect third order products Involving elastic variables, 
and fourth order products of small Quantities in^jj^qj^ 
following the expansion 0 = Oq + <9^* flatwi8e-edge- 

wise modal equation 


Ou 




4~ 


) Ik ^ 


(a) 


(4.65) 






is approximated to exclude second order products of 
elastic varieUoles, and third order products of small 
quantities, where it is understood that ajv, bj^, and 
apply for a zero torsion component 0^ in the flatwise- 
edgewlse mode , j . 

The result of these approxliaations is 


J K - ^ r f^K ^ ^ ^ K Pj ^ I a. ^ 

^ 


(4.66) 
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(4.67) 
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Coefficients (a^i*)o ^ ^®22*^o ^ (4.57) with 

substituted for $, Emd (a^i*)! to (a 22 *)i are equal to 
( da3^2*/d^) 0~^c> ®tc. Appendix 14.1 contains the explicit 
expressions for the latter coefficients. 

Recognizing that (Ax" ) " are the internal resisting loads 
due to vertical and inplane bending, the first term in 
(4.46) contributing tOy4^j]^ is seen to be a correction to 
the modal stiffness, bj]^ (bji in (4.7) ) , for differences 
in bending stif fnesscs''at the actual pitch, 0 , and the 
pitch, 9^*, used to find the normal modes. Equations (4.56) 
and (4.58) show explicitly the effects of pitch differences 
on the bending stiffness correction, Equation (4.59) 
shows that the second term in (4.46) similarly corrects for 
the uifference between centrifugal stiffenings in pitch. 


4.3 Addition of Structural Damping to the Modal Equation 

To represent the effects of structural damping in any mode, 
modal equation (4.7) is replaced by the equation 




(4.69) 


C 


■J 'C' 


= structural damping coefficient 


3o 




"j 


I ^ 


t h 


We show below that this form has the desired properties of 
structural damping. First it is emphasized that , w, v 
are elastic displacements and, consequently cure negligible 
in essential rigid body motions, e ''her pitching, f a. oping, 
or leading, with the result that is negligible in 
essential rigid body modes of bl::i( .notion and the structural 
d 2 unping correctly disappears for such motions. In essential 
elastic modes with weak rigid body participation 
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^ ^ 4 J 4^ * vv^ ^ v,^] 


(4.71) 

taking into account the expressions for w^ and vj^ in 
(3.6). In the free vibration of such an essential elastic 
mode, (4.69) reduces to the single degree of freedom 
equation 




V UJ Q. 


K I 


- O 


(4.-,2) 


in view of the properties of ajj^, (3.60) , and the damping 
term is the conventional viscous damping equivalent to 
structural damping, thereby justifying the term 2cc-i)^qjt 
in (4.69). 

A more tractable form for the damping term is one employing 
the critical damping ratio g* (twice the value of c/ccrit) , 


/^CRIT- 


(4.73) 


^r^lT ' ^ 


(4.74) 


(4.75) 


and (4.69) becomes 

^Jk ffc ^ V = S ,4.. 
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No provision is made in the program for input of different 
modal damping ratios, and only one value of g* may be 
loaded to represent all modes. 

With structural damping included eluent tj in (4.16) 
becomes 

tj - Qj - Y'jk ) 7/^ 

(4.77) 

and this is the current form in the program. 


Organization of Subsequent Chapters 


The subsequent text is directed to the finding of the S 
and T matrices in (4.20), and this section describes the 
organization of this work. 


We have already found contributions ajk# and to 
the elements Sjv and tj of the S cuid T matrices, tsj]^ 
is defined in \4.15) and tj is defined in (4.77) ). "'The 
remaining effort aims at forming Qj which will complete the 
sj}c and tj coefficients. 

The first step taken is the selection of the generalized 
coordinates appeeuring in the equation of motion in physical 
space, (4.1). These generalized coordinates aure the 
translations and rotations of rectangular axes defining the 
displacement of a point on the blade. A chapter is allocated 
to the definitions of rectangular axes. 


This is followed by a chapter yielding general expressions 
for the generalized forces in physical space, Fi, F 2 « and 
F 3 required for the formation of the generalized force in 
modal space, Qj ( (4.23) or (4.24) ). These general 
expressions for Fi, F 2 , and F 3 involve blade element loads, 
which subsequent chapters are dedicated to finding. 
Applications of rectangular axis transformations defined 
in Chapter (5) come into play here to calculate bending 
strains and bending reaction loads entering moment and 
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loading equilibrium equations. This ground is covered both 
to establish F? and F 3 and to emphasize that the scope 
of the analysis ^braces large pitch angle and twist 
rate terms without restrictions. 

Distinct chapters aure allocated to the formation of 
inertial, aerodynamic, Q^, and point loads generalized 
forces, ■' 

To derive Qj^, rectangular axis transformations are used 
again successively to derive blade accelerations and these 
yield inertial loads for substitution in the Fj^^, F 2 ^ and 
F 3 “ expressions. 

We then derive aerodynaunic foru-^s amd relative flow velocities 
on which these forces depend to form and this is 

by formation of the point loads generalized force 


followed 
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5. Rectangular Axes Detining Blade Displacement 

We define here the rectangular axes employed for the 
description of the displacement of a point on the rotor 
blade, including elastic deformation. Clear definitions 
of these axes are important because the axis transforma- 
tions, comprising linear displacements and rotations, are 
the generalized coordinates in the equations of motion in 
physical space. 

We restrict the text here to thv’ description of each 
individual axis which we define by means of a transforma- 
tion matrix giving the position and orientation of the 
axis in terms of a preceding axis. Successive multiplica- 
tions of transformation matrices to obtain formulas for 
displacement, acceleration, aerodynamic velocities, and 
force components occur in the text of chapters following 
this one. 

A complete description is given here of all coordinates 
needed in the blade equations of motion. We do not describe 
the support system generalized coordinates used for the 
rotor with grounded support, or for the rotor coupled to a 
rigid body in free flight which are described adequately 
in References ( 3 ) and ( ^ ), respectively. It is well 
to understand that these coordinates together with the 
blade generalized coordinates are the complete set of 
coordinates for the coupled problem, and that the present 
chapter is insufficient as a description of the complete 
set of generalized coordinates. 

Our blade coordinates resemble closely those used in the 
analysis of rotor blade response employing uncoupled 
assumed modes (Reference ( 1 ) ) . Similarities and 
differences are pointed out in the text. 

As a preliminary, we note that a compact presentation of 
definitions of rectangular axes is achieved with matrix 
notation. A columii vector of coordinates is expressijd 
as a capital letter with an overbar. For example, X = 

X, y, z. Trans forrxation matrices relating one set of 
coordinates to another are denoted with capital A and ^^6 
subscripted to show the rotation angle. A typical rotation 
transformation matrix is 
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0 


0 

/ 

0 


0 



(5.1) 


The symbol A_^ indicates that the rotation angle argument 
is . The relation = A_^ is extensively used in 
the text of subsequent chapters, and the subscript notation 
is valuEdile both for specificness and for easily deter- 
mining the inverses required in the derivation. An example 
of a transformation defining the relation between X 4 and 
Xr rectangular axes is 


^ ^ ^ 5 


(5.2) 


which is equivalent to 


Axis Definitions 

The ultimate reference is the stationary inertial rectan- 
gular coordinate axis, .'"ixed to the ground, whose vector 
of coordinates is Xj, Figure 4. Two different shaft 
oriented axes are used to treat respectively the rotor 
with grounded support, and rotor coupled to a rigid body in 
free flight. 
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Figure 4. Shaft-Oriented Axes for Analysis of Multi- 
Blade Rotor With Grounded Support. 
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1) Grounded Support Shaft Axes, X i 

Figure 4 illustrates the shaft axes Xi used for the analysis 
of the multi-blade rotor with groundea support. The X, 
system is reached from one translation and three successive 
rotations -jMg , Oq, defined by the following transforma- 

tions 



1 



I / 



(5.4) 


(5.5) 


(5.6) 


(5.7) 


The X, system is identical to the shaft axes of Reference 
( X ) . Because the hub is generally accelerated, Xj is not 
inertial. Angles same as those used 

in Reference ( 3 ) for the multi -blade rotor with grounded 
support. Xj xs the location of the hub measured from the 
Xj origin. Rotation transformations A_^, A_^g are 

listed in section 5.1, as are the other rotation trans- 
formations given below. 


2) Rigid Body in Free Flight Shaft Axes , 


Figure 5 illustrates shaft axes ^sed for the multi-blade 
rotor coupled to a rigid body in free flight. One trans- 
lation and three successive rotations yield 



(5.8) 

(5.9) 

(5.10) 

(5.11) 
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Figure 5. Shaft-Oriented Axes for Analysis of Multi- 
Blade Rotor Coupled to Rigid Body in 
Free Flight. 
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Figure 6 shows the relation between shaft axis and 
shaft axis The same relation applies betw^n X| and 

In algeSraic terras 

= (A-e' )i>f X 

Angles^', , y' defining the altitude of the rigid body are 
the sarae as those of Reference (2) . 

3) Rigid Blade Axes , X^^ to Xg 

Figure 7 shows these axes. They are described as rigid 
blade axes because they would yield the displaceraents of 
a point on a rigid blade involving azirauth, , hinge 
offset, e, radial position, r, and rigid blade displaceraents 
6, , 0, only. These axes are identical to axes with the 

sarae ncuaes used for the rotor response analysis eraploying 
uncoupled assuraed raodes of Reference (1) . They super- 
irapose displaceraents on those due to shaft displaceraents 
and are defined by 


X, - 

hr 

(5.14) 



(5.15) 

^ b - 


(5.16) 

X4 


(5.17) 

X = 


(5.18) 

- 

Ae 

(5.19) 

- r 

h 


(5.20) 


0, 0 
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Figure 6. Relationship Between Shaft Axes Used for 
Grounded Support and Rigid Body in Free 
Flight. 
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4) Blade Elastic DiaplacamentS / X 7 to X^q 

Figure 8 shows the axes superposing elastic displacements 
on all the displacements described above. Elast^ dis- 
placenents are u*, v^, Wg or alternatives *• 

u*, V, w ,amd elastic twist B . Displacements v- and Wg, 
or V and w, cause angular deformations of the blade A and 
^2 which are expressed in terms of Vg and Wg or v and w 
in section 5.2. In addition a displacement » 0, v», 

Wq accompanies the elastic rotation of a bent elastic 
axis, as shown in section 5.3. These W. dlsplacments 
are shown in the text of subsequent chaptexs to be equal 
in Importance to terras due to w, v, 0^ alone and are 
justifiably included. 


The first difference between the axes usad here and those 
used for the euialysis of Reference (l)_employing un- 
coupled assumed modes is that the new X-j axis differs from 
the uncoupled modes_Xio axis of Reference (1) by a trans- 
lation Wg fr<xn the Xg axis. The axes are defined by 

- X-J W,' (5.21) 

^7 - ^ \ (5.22) 



A Y 
A| cj 

Y f W7 


with 






w 


(5.23) 

(5.24) 

(5.25 



0 , 


w 


(5.26) 


Coordinates W are preferred to for the final equations. 
The latter, it is recalled, were the elastic displacements 
of Reference (1) . Defining 


l/V 



V ^ u/ 


(5.27) 
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Figure 8. Blade Elastic Displacements for Mult^-Blade 
Rotor Analysis. 
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where v and w are illustrated in Figure 8 and letting 


X 

' \aJ i R 

(5.28) 



(5.29) 

/X ^ 

- r 0 0 

“ / ' 

(5.3C/ 

we 

find from (5.18) that 


W 


(5.31) 



(5.32) 


and (5.21) may be replaced by 

■- X 7 ^ A W (5.33) 


5.1 List of Rotation Hatrices 

Listed below are rotation matrices required in equations 
(5.4) to (5.33), 

(5.34) 

^ f'S ^ 

- O 

o 0 I 
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( 5 . 41 ) 


( 5 . 42 ) 


( 5 . 43 ) 


( 5 . 44 ) 


( 5 . 45 ) 


( 5 . 46 ) 
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(5.47) 


Matrices A>2» AAi» h0^ result fr<wi a small angle approxima- 
tions to angles \2> (Figure (8 ) ) , such as was made 
in Reference (1 ) . 


5 . 2 Derivation of Angles \ and A ^ 


Figure 9 illustrates angles emd A 2 « It is scien that 


■^3 ' /d2-v ) 

(5.48) 

' ■ ^^^7 / ^ Ci 

(5.49) 

From (5.18), (5.33) we find 


^ - '^(9 ^ ^ W ) ^ 

(5.50) 

' ^(9 ^ y ^ W 

(5.51) 

^7 = (x --CC/ - 7^ ^ 

(5.52) 

" (x-'W~R ) 

(5.53) 

Writing (5.53) out we find 


^7 - 'X-f'4u,^ 

(5.54; 
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= ldy,/dx, 

X, = (dz, /dx, 

Figure 9. Derivation of Angles Ax andA 2 * 



Figure 10. Elementary Twisting Couple Acting on a 
Blade Element. 
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(5.55) 

(5.56) 


The changes in xn, yn and z,, for a variation such that 
y and z lie on the elastic axis, are 



= ( 1 4 u ^ ; d)z. 

(5.57) 

1 1 

(Lty + 5^<1W 

(5.58) 


- dv ^ cTw 

(5.59) 


It follows from (5.48) and (5.57) to (5.59) that 




/ ' iX 7 ^ '• 




A/ 




(5.60) 

Prom (5.22) 



liv-, 

" ^8 

(5.61) 



(5.62) 

( i ^ ■ CL 



^ ^ ^<9 ^ 

(5.63) 
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In terms of the and w^ coordinates of Reference 

these results are 

( 1 ) 

a 

/ 

-• ^ ^'e + ' 

(5.64) 

A , 

y \/^ f u/(o 

(5.65) 


which agree with (1.23) and (1.22) of Reference (1) 


5 . 3 Translation Wg Accompanying Elastic Rotation 0 ^ 

In this sec .ion we derive the translation of a blade section 
located at X 5 = x accompanying elastic rotations 0 ^ occurring 
at stations radially inboard of this section. The trans- 
lation of the outboard section is represented as a displace- 
ment W"" =0, vg. Wo of the origin of the Xin axis from the 
X 9 «ii, <Fi 9 url 8 )° 

To calculate vg and wg, imagine tlie blade to be loaded by 
a twisting couple acting on a small length of blade. Figure 
10. If the inboard section of the blade element is fixed 
and the outboard section rotates by dd^^ the entire blade 
between the tip of the blade and the outboard section of the 
blade element twists through d^. Rotation occurs about a 
straight line passing through the local elastic center of 
the element to which the elementary couple is applied. 
Superposition of such elementary couples yields arbitrary 
loadings emd corresponding ti^ist distributions needed to 
calculate elastic center displacements accompanying rotations. 

Figure 11 shows the relatl/e positions of a straight line 
xq directed along the elastic axis at x, at which the 
eleaentary couple is applied, and the elastic center Q at 
section x. P is the point of Intersection of this Xg axis 
with the y - z plane at x5 » x. To a high order of approxi- 
mation the rotation arms about this elastic axis are ^y* and 
hz*, shown in Figure 11. Neglecting the x- displaceaent due 
to rotation, the displacements of the elastic center due to 
the elementary couples at x are (Figure 12) 


PAGE 




Figure 11. Geometry used to Calculate Elastic Center 
Displacement Due to Elastic dotation. 
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Figure 12. Geometry Used to Calculate Elastic Center 
Displacement Due to Elastic Rotation - 
Concluded . 
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CX 






(5.66) 

(5.67) 

(5.68) 


The coordinates of P are 


V i- V ( X - >c J 
A ^ /\ \ 

Z = w -t- W \ Y - -X ) 

and 

-V-tj rV-»/-v 

^ A /V ^/ -^ \ 

Az - w- z -w -\A/-w(x:-i>cj 


(5.69) 

(5.70) 

(5.71) 

(5.72) 


Hence 

S. -- 



(5.67) and (5.68) become 

K/ - C- V (v- V )1 6lSr 


(5.73) 

(5.74) 


Superposing all the elementary couples from 0 to x to 
generate the actual twist, we find the displacement of the 
elastic center at x due to blade rotation to be 
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y 13. 

A V ~ J - wCv j - w Cr ^ 


The corresponding displacements in the - Zg plane at x 
are obtained from a consideration ofthe^ transformations 
relating free vectors in the Xg and X axis systems. Denoting 


(J V , 


(5.77) 

(5.78) 


we find from (5.18), (5.21) to (5.23), and (5.33) 

(y X ^ A ^ X cj 


or using the reversal law :or reciprocated products 


cH 


M 


^ 'A, ^ 'A 9, 




-- ( 


-A 


from which it follows 


- V 

*^(9 
- 5e 


^0 


(5.80) 

(5.81) 
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(5.82) 

0 

Zl ^ 

4 0^ A V 
6 

(5.83) 


^0 A. 

(5.84) 


These are the displacements of the elastic center at x due 
po elastic rotation of the blade. Hence, the translation 
of the origin of the axis from the origin of the X 9 
eixis is 


(5.85) 


with (5.83) and (5.84) yielding V 9 and W 9 . 
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6. Moment Equillbrlmn Equationa and Expressions for 
Generalized Forces Fir 

To form the generalized forces In physical space Fi, F2/ 
and F3 required for Insertion In the expression for the 
generalized force in the modal coordinates, Qj, we first 
derive internal resisting moments and equilibrate these 
to the external moments to form moment equillbrixim equations. 
Comparison of derivatives of these moment equilibrium 
equations with the general equation, ( 4 . 1 ), yields F^, F2, 

F 3 * 

To calculate the internal resisting moments we first 
calculate the longitudinal strain of a fiber of beam amd 
then sum the moment effects of the strain in. the section. 

The next section provides a very brief description of 
the derivation of strain, and this Is followed by a section 
yielding the Internal moments. Expressions for 7 ^, F2, and 
F3 2u:e then calculated and this is followed by a section 
on expressions for moments ana their derivatives occurring 
either in the F expressions or required for display of the 
system response. 


Strain Expressions 


To calculate the strain component perpendicular to a 
norm,!,, section of the beam, we consider two adjacent 
secti'>AS, separated from each other by a distance dx© ■ 
dr in the unstrained state, and calculate the ext Vision 
of a fiber extending between the faces and passing through 
points at yio ^10 in the faces. The usual ass'imptlon 
that plane sections of the beam remain plane after elastic 
deformation is contained implicitly in the geometry of 
deforma' ion assumed, and described in the previous chapter. 
The expression for strain is 


e - 


d/_y, 

cOs^ 

^'Yo 

<^So 

^0 


d^o 


( 6 . 1 ) 


( 6 . 2 ) 
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where ds » fiber length after deformation and ds^ » fiber 
length before deformation. If we express strain^in 
terms of components of position in the Xe (that is, X) 
system, then we must evaluate 



(6.5) 


Repeated application of the transformation's of the previous 
chapter yields 


X' 








( 6 . 6 > 


This expresses_the_X displaccmen'. j ir terms of elastic 
displacements W , and elastic angles eind A 2 (related 
to Wfl by (5. 64) ®^uld^ (5. 65 ) ), and section or rdinates R 
and Xj^Q. 

By asssiming all elastic dispxacements to N absent in (6.6) 
we obtain the cc^ordinates of the point before elastic 
deformation. 
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0 - ^ I 

^ ^ lo 

After elastic deformation, we find 

u* - A, 2 

-6^ 

z -- 2 ^ . ®e Ce •Se ^ “^e^e ^ ^6 


To derive (6.8), 2nd and higher order products of elastic 
variables were neglected. Carrying out the differentia- 
tions of (6.5), we obtain 

' .- u' '+ T- Aj'i- We'e-' (- V(, (0')^)^ ,o 

e\\^e(9'f) 2,0 

+ ©'e/Cz,o ^ A,o ) I6-9) 

Substitution of the expressions for and A 2 ( (5.64) and 
(5.65) ) in (6.9) yields 






- \/(o - G - 3yv/^ G ' ^ ^ 

- Wf7 ' f (9 y. 3. ^ 

9 0/ 


( 6 . 10 ) 


C 
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6.2 Internal Resiating Forces and Moments 

In this section, we calculate the loads of reaction induced 
by elastic deformation, employing strain expression (6.10). 
Figure 13 illustrates the positive conventions of these 
loads. The positive direction of Myj^n is along yj^Q and 
is not opposite to yxo ss was the positive moment employed 
for the determination of the normal modes. 


Referred to the x, y, z axes, the strained fiber has 
direction cosines x'/s' , y'/s' , and z'/s'. The corres- 
ponding tensile reaction on the face of the beam has 
ctxaponents (x'/s') E dA, (y'/s') E dA, and (z'/s') 

E dA along the x, y, z axes. The corresponding elementary 
reactions referred to the Xxq system are 



A 


-Q ^ ^ ^ ‘ 

-Q 


" A T * A . 




-A - 


-9 


- s 


e 


> v>/ 




( 6 . 12 ) 

with 8 being the total angle, defined in (6.34). From 
(6.11) we obtain 


d2 ^ cC'A 

(6.13) 

^ ^ ^ ^ ^ ^ ^ ^ ^9^ ) 

(6.14) 

Ct' Fz /o - ^ \ ^ Cq Z ') 



77 

PAGE 





I 


i 


Sikorsky PIrcraft 


u 

DIVISION UNfTSO AACilArT CONSOAATK3N 


REPORT NO SER -50912 


In these expressions second order products of small 
quantities involving , y' » z', A 2 and are neglected 
to obtain (6.13) and third order products are neglected 
to obtain (6.14). 


The section moments are 


M 


to 




dl'Tc 


/O 




•a 


I O 


^ b /o ^ Z,, ^^,0 


(6.15) 


(6.16) 






O/o ^ 


Moment Mxxo is evaluated at the elastic center with the 
result that no contributions to moment come from shearing 
stresses, because the elastic center is defined to be a 
position where such torque is zero. It is not necessary 
to approximate dPxio# (6.13), to a higher order than 
indicated because Myxo Mzxo neglect second order 
products of elastic variables, and third order products 
of small quantities. The approximations to dPyio and dFzxo 
are high enough to permit the consistent determination of 
Mxxq to an order which neglects third order products of 
elastic variables, and fourth order small quantities in 
the torsion moment equilibrium equation. 

Substitution of €, (6.10), in (6.13) to (6.16) yields 

M X/o - G J r ( (9 f 
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I 





' 

" ^ 4 - i- A V C S’ 

^e(6 'A 

Hz . 

-e^ Px,, -C-, tx. 


- 

'"a 

Cl 




with 



^ r € EcM 

’'p\ 





tr, : 



■- 

^ 'j ,o A 



h 1^,0 


(6.18) 

(6.19) 

( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 

(6.23) 

(6.24) 


The above results were obtained for sections with structural 
symmetry about the yio axis. 

The underlined terms eure additional to these present in 
equations (1-51) to (1-53) of Reference (1) , and represent 
the added effects of large twist rate and non-linear twist, 
without restrictions. 
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To represent a Sikorsky type oounter%^ight we add to the 
right-side of (6.19) the term 



c w 




(6.25) 


where <^eAcw distance between the chordwise position 

of the centroid of the counterweight and the elastic 
centroid at yio » e^^ , and r^^^ is the inner radius of the 
counterweight. The loading (PX]^q)cw the inertial load 
of the counterweight per unit span. 



(6.26) 


In (6.26) m_y is the mass of the oounterweight, m is the 
mass per unit span of a blade section, including the 
counterweignt and px$g inertial load per unit 

span due to mass m. 


With the counterweight term the expression for Mzxq is 


/V 


7c 


^ fA 


f ^ '/• 


/o 




j; 

' o 


r"K 


C.<jJ 


Vn 


^ H 


+ ( 


/; 


C UJ 




C< 5 ') 

(6.27) 


TO be ^d>le to identify later the generalized force F 2 and 
F 3 in vector F of (4.1), vdiose last two rows express^ 
vertical and inplane equilibrium of_loads, we transform 
structural reaction moments to the Xc system and using 
(5.32) we replace w« and Ve by elastic displacements w emd 

V. 

By means of (5.32) we can show 
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^(9 ^ ^ \/C= a 


II 


^0 ■* '' 


(O'J’ Vg- © ". 3 




Vv^p f 


Substitution of (6.28) and (6.29) in (6.18) and (6, 
yields 

-- -e I. 0,, V " * tx , s„ V " 


io ^ s ~ - z e 






■^a' x-,o 


Dcuj 


Moments are related by 

in,- = Ag/V,A^, 


with 






W 


■> ^ B 


k 0.-A 


S 


C 


G 


e 


and 


© A G 4 &, 




G ^ 9 
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(6.30) 

(6.31) 

(6.32) 

(6.33) 


(6.34) 
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Expansion of (6.32) yields 




(6.35) 





(-M, \ 

J ( o J 


Q-^^z 


(6.36) 


To obtain (6.35) and (6.36) we neglected second order 
products of elastic variables and third order products of 
small quantities. 

Substitution of (6.17), (6.30) and (6.31) In (6.35) and 
(6.36) yields the desired vertical and Inplane equations 
for structural reaction moments in the X 5 system In terms 
of displacements w and v. Assembling these equations with 
torsion equation (6.17), we have 





(6.37) 





A 


A 





(6.38) 
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(6.39) 


Ct J " + 


(F I » + L- r c * \ " 


The corresponding loading equilibrium equations are obtained 
by differentiating (6.37) to (6.39). 


(6.40) 


-[C7&;Ue'*0e')<A' o 




V 


(6.41) 

II 


^0 


I ^ 


^ A <? ,, r 

''c^J 


p>‘. 




/o 




tri 


G 

1 1 1 

i)j F /1 


'Cu. 


3, 

V6.42) 


O 




^ A -^,0 ^ cw 


/ 


o 


Cj£J:^ i ) j ~ ' O 




I ^ 
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6.3 


Generalized Forces Pi and 

Examination of the rows of the general equation (4.1), which 
are 





3. 

^ (S 



r F, (6.43) 

(6.44) 

3.. y // 

+ ' 6 ) ^ _ 

(6.45) 


indicates that subtraction of (6.40) from (6.43), (6.41) 
from (6.44), end (6.42) from (6.45), will yield Pj^, P2 
and P3. Doing this and employing the properties wj^" « w" 
and vi" ■ V", we obtain 




|0 


) 6, 


)C, &e 

Jio ^ 


t- 


r « 

)e, 


/■ ' 


Kfc F.,J 




(6.46) 
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- ^ f (-Tv/. ) r>\ VA/ j 

■'\w 


ir>^ 

io 


C w 




J 

(6.47) 


h z. 5 -^ (- T v/ ) - r>> ^ V, ^ 

^ y '-pP^ i 


II 


C SA,^ 


(6.48) 


We express Mx£q, Mye* > Mz^' in terms of local torsion (qx5, 
qys, qz5) emd airecz (PX5, pye, PZ5) loadings per unit 
span to reduce further these derivatives to useful forms. 

In the next section, consideration of equilibrium of a 
beam element Is shown to give 


(6.49, 



1 

- V 

^ ^ f vv* ^ 



- w 


i 






-t V 
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where 


^ ' S V/s / 


(6.50) 




r'p 3 = 


r 


z 


J 


Vz 


A/ 


With 


n 


lo 


A , l\ 

"e 


-a '’-X. 


i -e b 


(6.51) 


and (6.12) for A-Aj^ A-A2 A-0 , %»e obtain 


(A 


Ac ■- W ' lA. ^ W>1 


^.5 


c 


(6.52) 


Differentiation of (6.52) yields 




fw"M, <«-53> 

s 


and substitution of (6.49) in (6.53) yields 




-V 


''“is ■ '''^ ‘'z. > 


V " M 




■\ Vn/ 


" H 


^ 5 

(6.54) 
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Substitution of (6.54) in the expression for P., (6.46), 
and use of the expressions for (6.38), ana 

(6.39) yields ^ 


f 




+ yv\ 




J \o 




+ lOA \ / 


+ (e'V' 

V/ eg ^ 



(6.55) 


Using the roL .its of the next section that 


M 


- -‘1z 



in the expressions for F,, (6.47), and F3, (6.48), 
obtain 


(6.56) 


(6.57) 

we 
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. + (-Tw/)' ^ 

- P p ^ £ ) 

C'c ^ irv^ -I 

(6.5 

Pj ^ -‘U; vCVF^,_)'v ^ C-Tv:)'-^n5i^v, r 

- [ce(-eAFx,.+ y^^*" Py” .iq^w. .g ^ )1 " 

(6.5 

Because we desire to express F^f F 2 and Fo In terms of 
X 5 syston loadings, we ultimately must substitute 

^ ^ 5 ^ S (6.6 




t- V w 


in (6.55), (6.58) and (6.59). Ma)cing these substitu- 
tions as convenience indicates, expazMiing 


G ~ G ~i-(y e 


and deleting negligible terms, as appropriate to torsion, 
flatwise, euid edgewise equations, we obtain 


•»«, I Fv n 
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+ ^'55- 

t W ,- 


-(( T-Fjry'jK^ fcj,') 

-^(g'k * Fx,„V 


*-Ce -^A 

Afo 

w 




- ^0® A 

Fx 

A, a 

I) 

V 

m 





D 

Px,, 


i)4 



CUO’ 


rvi 



r 

w ^ 

*~®c w 

po 

* y* 

^ / O 

I^C VA/ 
m 

-d)^: 

F fc p ( V 

M 

" c 

V/ j£) 

) (-^ 

(\ 



(^e 


i\ 


f P ^^cw' 

J ' X’j v>v 


dn 


o 


^\A/ 






V' S 


Z6 


(6.63) 
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V 4 ^*'a,„, rVy< c{ ^ )' 


( I 


'A Ov 

' ^ vy J I V s 

^'c ^A/ 


(6.64) 


^ ^'35- ^-f'T V//) • +l^W'^ 

- J-' &(-(’>, J f) 

r ' yv\ / 


■ fee 5a (- 


C^ 




A 


r "* ^ 

X 


// 


ew , 


nFN 


(6.65) 


The aim of the subsequent chapters is the derivation of 
the elements on the right-sides of (6.63) to (6.65) in 
terns of the coordinates of the problem, and the sub- 
stitution of Fxf Fa. and F 3 in the modal force expressions 
(Qj'^^» Qj'* » (4.26) to (4.28) ) to calculate Qj. 

The last term in the torsion iKxaent F^, (6.63), is the 
torque from bending forces (compare p.430ff« Reference (7) ). 
We verify that for » 0 (linear twist) and negligible 
the term is 

(ei,-ti.J(v s^e)- 

(ti (vp " wp"- a €* ‘wp " v/p '4 3 Q 'Vp 'Vp' 3 

( 6 . 66 ) 

which checlcs with the equivalent term in equation (1-70) 
of Reference (1) . 
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6.4 Beyn Elyayit Equilibrium Equations Used in the 
Multl-Biade Rotor Analysis 

In this section we state briefly the results of considering 
the equilibrium of an element of beam^ to justify express** 
ions for moment derivatives us*‘d in the last section, and 
to derive moment expressions used ultimately for display. 

Figure 14 shows the forces and moments acting on an element 
of beam. All forces and moments are assumed to act at the 
elastic center to be consistent with the derivation of 
structural reaction moments, and eire resolved to the 
system. The results of force and moment equilibrium 
considerations are 




1 1 

0 




( 6 . 

67) 




0 







V 

■= 

c 













( 6 . 

68 ) 



X 5 T V 


, - w 

0 



^'^ 5 - 



'fxs 

' 

— 

c 





1 

^5 ' 



- 

0 



Differentiation of 

the 

last two 

equations in 

( 6 . 68 ) yields 


- 


■ 

' - 



( 6 . 

69) 


- 

■lx,- 

4(V 

' + 



( 6 . 

70) 


which are (6.56) and (6.57) used in the last section to 
obtain Fj and F 3 . 
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Integration of (6.67) with respect to v and insertion 
of boundary conditions that tip loads are zero yields 



r- ^ 

which is (6.50), employed in the moment derivative 
expressions of the last section in the derivation of the 
F forces. 

To obtain moment expressions Mx^, and MZe> we note 

the following identities. ^ 

V'(h) --/£[(v(e)-v(x))F^^(f)] (6.7J) 

+ [vCf ) (fi) 

- ('w(f') 

vy'(£) Fx3 if) (>f)) ), 

ff w(^ ) - W ( X ) ) P X 5 U) 

~v'(f}F,Jf) ---X [(v(^)-v(x^)fx,(^) 
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(6.73) 


In (6.73) C| and Cy are arbitrary constants. Identities 

(6.73) may to verified by differentiation. 

Integration of (6.68), substitution of (6.72) and (6.73) 
in the resulting integrals and employment of the condition 
that tip shears and moments vanish, yields 


y 5 -- j* 




J5" 


’ J T'~(^)-w(x)] )cf/ 


(6.74) 


[w(^3 - o-'i' 




(6.75) 
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-£u^) -(^1] 

^ f 7 ^' 

(6.76) 

The expressions for local blade oriented moments derive 
from (6.51) with (6.74) to (6.76) supplying the X 5 moments. 



^ X ^ t- , V- Vt' ' M2 ^ 

5 j b 5 

(6.77) 


-Aj rtAj Cg f 5g 

(6.78) 

^/u ~ 

-A, >-C^ 

(6.79) 


To obtain (6.78) and (6.79) we neglected products of 
elastic variables, -Ai 0 qMx 5 andA 2 ^e^ 5 * 

Moments (6.77) to (6.79) employing (6.74) to (6.76) are 
termed the external moments inasmuch as they derive from 
external sources. Moments (6.17), (6.30) and (6.31) are 
termed internal moments because they are the equilibrating 
Internal resisting moments. 
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7. ACCELERATIONS 

We derive in this chapter the acceleration of a point on the 
blade, needed in the expression for the inertia contribution 
to the generalized force, Q^. 

In our ^proach, the acceleration components are referred 
to the Xe system, and these accelerations are required to 
involve ue linear and angular velocities and accelerations 
of the hub in the shaft oriented Xi system. The pareuneters 
of the Xi system motion, Vqx And (Figure 15) and their 
time derivatives eure assumed known. These pareuneters derive 
from the solution to the coupled support or coupled rigid 
body system of equations, obtained according to the follow- 
ing procedure. 

The response of the support (or coupled rigid body) derives 
from the input of rotor forces and mo’^ents to support ^ 
at the hub. Then the support responses, v^ 
are returned to the rotor blade equations o 
rotor forces euid moments are calculated. These new forces 
are fed back again to the support or rigid body, and the 
process continues with time to obtain the response of the 
combined system. Figure 16 shows the flow of calculation 
coupling the rotor and support modules. 

In view of this procedure, our plan is to express the blade 
acceleration in terms of the motion of the hub in the X^^ 
system, as well as other motions, and to resolve the accelera- 
tion components to the X^ system. 

Following this, we elaborate four different versions of the 
acceleration with a mind to easing and systematizing the 
subsequent derivations of inertia loads and generalized 
force Q^, and to effect a separation of modal accelerations 
q from q and q in the modal equation. 

Including the acceleration of gravity, the acceleration of 
a particle on the blade is 


^ 1 ' ^ol ^1 
I motion and new 


^ c7 

cCtr ,'i 


(7.1) 


where the velocity of the particle may be expressed 
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Figure 15. Rectangular Axes Used in the Evaluation 
of Acceleration. 
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ROTOR HEAD 
FORCES AT t 



Figure 16. Flow of Calculation Illustrating Rotor/ 
Support or Rotor/Rigid Body Coupling. 
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(7.2) 



+ 


di t- 


(7.3) 


In (7.2) and (7.3), and Vqs are linear velocities of 
the origins 0 ^ and O 5 of the_Xx and X 5 auces, is the dis- 
pla cement vector and P is the displacement vector 

OcP defining the position of a point on the blade, Figure 
15. 

Substitute (7.2) and (7.3) in (7.1). We obtain 


^ ' cX t 0 i 

' 2 — 
£ ^ 



dt^ 

(7.4) 

— _ d< 

' oit 


(7.5) 




■ rft' ^ 

dit 

(7.6) 


The grouping in (7.4) eases_the evaluation of Integrals 
Qj by distinguishing te^s a^^ g, independent of blade 
radius, from the term d^P/dt^, which depends on radius 
through radial variations of pitch and elastic displace- 
ments. 

— — 2 — 2 

We proceed to evaluate now a.c, g, and d P/dt , which 
compose the acceleration a, T/<4). 

a) Term i ps 

To obtain the acceleration of the origin of the X 5 system, 
Sq 5 , we express 


( 
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(7.7) 


CL 


O . 


and using (7.5), we find 

r ^ A 

(£i . 


K-3.5., ^ J, ^ 


^ 4>' 

<ict 




c(^Ar, 


(7.8) 


Employ the properties 

Ctt, .- - 

- I \/ 1 

2t 


- ou, yi L f ^ -j 

^ 'jt 


i' -- ’‘J, ' ^ ‘'J, 


where 


C*_> r L , o.>v f K z 

f / ^ , J I ^ I I I 


(7.9) 


(7.10) 


Define 


We find 


,/S’ 

(s-c 


(_ P( . . d 


>at 


-t K 4 


a 






- & -v), 


+- UJ ^ xv^ 


cQt 


S' 


Simileurly, from (7.3) we obtain 

'^0^ -- Cb, + f£|L), + C^, X, 


/O 


(7.11) 


(7.12) 


(7.13) 


where 


VC 


0, = S ^ s, ^ A;, 


w’, 


0. 


(7.14) 
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P -- w ^J, ^ 


(7.15) 


Define 




(7.16) 


We obtain from (7.12) and (7.13) 




“ 1 - 

(v^ ) 
! 5 i 

^^<55- ) *3 

0 

X-C 

- (v^.) 

- a. 

‘x, 


r° 5 ')l, 

= (v.,^ 

t t-Jx, C«o,-)u -‘^•>■^,('''05^1 

Wo^)xi ‘ 
/ . 

- ^<^x. 


ftz, - 

, /’^i 




,' ‘"’x, pz, 

/-X, 


: 

, ^1 


. ■ ^3 . , 


, (7.17) 


Substitute (7.18) in (7.17). We obtain 


) - V/) ^ 


t ‘ 


'3.' ~ 


(7.18) 


(7.19) 


r Lj - Lu^ p 


t -f -h / i \ 

■3: -u-^^ ^ 

'7_ i ^ ■ 


V ^'-3, '*'^3 I ' ‘-^<1 - ,/l, _ ^ <--Z z , ^ 
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(CL, 






I ^ X, VV, + ) 


To find /’, we apply successively (5.14) to (5.17) amd 
obtain 


/ 




(7.20) 


Since /* * Xj, for Xg = 0, we get from (7.20) 
Since B » constant 


(7.21) 


- /'V ^ (7.22) 

/T.-Ay'- (7.23) 

We obtain the following results 




/ 

r a7 

V ^ ^ ) 

0 ^Y , 


, /> 

/O 





Az 

/ 

. 

1 1 

^ y 

r^, 


,A 

/ 

- - <:• Y 

r. 2 


o 


,2 


^ 5 , M.' 'j 


(7.24) 


Equations (7.19) and (7.24) complete the formulas defining 
aiQ 5 in ^e system. To resolve the components of Sq 5 
to the X 5 system, we define 






(7.25) 
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and apply the transformation 




. \ 

^ (7.26) 

/ 








Using the repeated suffix convention for summations, we 
write this 



where the repeated suffix j indicates summation from 
j s 1 to 3, and where 


J J ^ '/9 ^ 7 ' (7.28) 

- >T3 ^ ^ S 


(7.29) 


; 0 


'5^ (. / 5^5% , • >'3 f 5- i'l,, • )<3 ^5 

, / \ /- \ /-V 


^ J 


t : /, S,3 


The symbol for the direction cosine, Sj,j in (7.28) is not 
to be confused with the generalized mass in normal 
coordinates , (3.58). 


Equation (7.27) defines the components of acceleration 
of 0^ in the directions, in terms of the motions Vp^ 
and uPl of the Xi system, and the generalized coordinates 
/^, S, and azimuth >/, and is thus the desired form for 

*o5* 
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b) Tarn g 

Gravity acceleration components are derived below for each 
of the two shaft-oriented axes. 


Transformation (5.4)_ to (5.6) yield the relations between 
free vectors and Xj in terms of the angles of the 
grounded support 


X 


I 




I 



(7.31) 


Identify the free vectors as 





(7.32) 









We 

obtain 





- ‘i 



o 


(7.33) 



'5 



To obtain the Xi gravity acceleration components in terms 
of the angles oi the rigid body in free flight, we use (5.12) 


Replacing by (5.8) to (5.10) this becomes 

/ 

-Q ]-jY ■ q' ^ 

and using (5.13) %ie get 


(7.35) 




- ^ 





(7.36) 
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■' - 9 ) ^'4 

Identify 

^ -- K , > *3 ^ 

We obtain the systoi components of acceleration fzxxn 

‘Sf.f - « 5x, = ‘'■'j -ij „. 3 „ 

3 b I ‘'■•'j 

‘yx, 

with S 4 j defined by (7.29), and the repeated suffix j 
indicating summation from j » I to 3. 

Equation (7.39) defines the components of acceleration 
in the Xc system in terms of tte X, system gravity 
acceleration components. The X| system gravity 
acceleration components derive from (7.33) for the grounded 
support and (7.37) for the rigid body in free flight. 

c) Term d^/dt ^ 

2— 2 

We derive here the acceleration contribution d V/dtr 
in (7.4), arising from the displacement, P, of a blade 
particle from the origin O 5 of the X 5 axes. 

Defining 


and differentiating this twice, we obtain 


U 

eijgmmT ( 

Pi 
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This yields 


^ i (^9 )rt ^ 


o 

c 

3 


(7.37) 

- 




I - ^ , 5 (7.38) 
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f? 


d "'T 

Ci' (- 


X 


• J U 7'P 

L - ^ ^ 

- ' 

_a ( 

‘ 2 ^ r a 7 X . 


^ Jy 


a 






/ 


Ot *• ut 

h I ^ o' ^ -r^ . f ri ^ 


^ J y' Cl' ^ i>' 2 " 

^ ^ f'»L 


fv ^ ^~ 


ai^ 


rci 


cx-i 


. y Jf )r, ^ ( \ 1^ 

rt' -^ / X, I -Tirr-?./ :?o I 


(7.41) 


To reduce (7.41) to a working form,.^ evaluate first 
dXc/dt, dJs/dt, ^ 5 /dt, d^ /dt^, d^-^dt^, o^ 5 /dt^ and 
this is folloirad by the expressing of^ in terms of elastic 
displaconents and pitch angles. 


Employ the properties 

— - cwT's- 


- tx- 


a 


v j,- , 




r X ^ 


r 

(7.42) 


Define 


Use the identities 

■>' -^5 


(7-. 43) 

c , . c 

J,- '' ‘■y - 

K^. y.Zl- 

"J y 

(7.44) 

‘-^/ 

/ 

\ 

It 

r b- - 

kj ^ 

' 


^ > — ^ 

C ' t ■ ~J v' 

y 

J - «- XT 

f<iT ^ 

- 0 


We obtain 

ckr - 

'• 0 

'■'TIT" - 

J 

dU ' ^ 

*' 3, ‘^^s- ' 

( C f ' 

- ^L, - 

J :? 

UJ .. ^ 


(7.45) 
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Na darlva now tha tan d^^/dt^ 



This requires the evaluation first of 
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The first tsjrrn in (7.46) beconss 


c^i: 


Kc 


c . -h ^ X 


r ^ ^ J 3“ > Tj -5 


, o -» Ji' ^'^5 - S' 


5 * 6 


(7.49) 


(7.50) 


after use of (7.44). 

The second term in (7.46) is 


,ic^ 


CO 


5' 








(7.51) 


after substitution of (7.43) and (7.45), and reduces to 

^ ^-2 ^ 

(7.52) 

Substitution of ^7. 501 and (]^52) in (7.46) and similar 
treatments of d^Jg/dt^ and d^^/dt^ yields 


cr Y 


d. 

ci t 




o^< 


‘1' 


*- CO 




CO, 




109 

PAGE 


.. i . 


- ^9 9i?V D 



1 


r 


1 


Sikorsky Pircraft 


u 

AmomAPr 


REPORT NO. SER-50912 




;~2_ 




(c . c 

^ jV-co^^ 3 ^ c * c., ) 




d ’"*/<i _ 

.-• 


cK- 

\ 



-•■‘^•x\- - <-' ■Jj,- '' '^’) 


Define 


d ^P 


- — -- ^ ‘i- 




Substitute (7.45) and (7.53) in the expression for 
d^P/dt^, (7.41). We obtain 






('Uj, 




^ )T-^^-^ 


I I ‘ 40 Ary n 



REPORT NO. SER-S 0912 


Sikorsky Pircraft 


bivtatON cm uNtno 



vnoM 


f “ X,' (^ 2 . ^3 y ) ^2 5 " 


•t- ( - ) Fz 


These acceleration components are reduced further to 
useful forms expressing ^5 and P in terms of the 
generalized coordinates and the motion of the Xi axis. 

The angular velocity of the syst«a is 


- - r, " <i — r^, 1/.36) 

^ 5 ' ■" ^ ^ ^ ' J r 

which, with all contributions resolved to the Xe directions, 
is 


/aJ . z 

'/3 - cS - ^ ^ -’<1 ^ -p 

X 0 4- 

G 



cJ . 

c 

■/^ 



^ 2 / ^ 

'S 

0 

(7.57) 


Taking account of the definition of the direction oosines 
aj,j in (7.29), the expression for in (5.43), and 
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defining 




equation (7.57) reduces to 


( ' / ,-7 

^s>' ■ - 7 ^ 

<^Z^- - ^'f3 ^ 

vdiere the repeated suffix j indicates summation from 
j = 1 to 3. 

The angular accelerations in (7.55) are obtained by 
differentiating (7.59) with respect to time. 


- o./ u; ^ ^ ^/3 S ^ 

j j J ^ J j , ' 


(<) ■ (0 " 


z 1. W tp 5,-3 .- Cp 5' 

J yj \J 

\ppendlx (14.2) defines coefficients j , which are the 
derivatives of ais defined iu (7.29). The angular 
accelerations are 

o) . , 

^ U)^ . CO . . - / ^ 


'^C - , ‘-'s, ; 

Obtained from (7.58), with « 0. 


/, S, i 


Equations (7.59) and (7.60) fnvoive the motion of the Xi 
axis, ^yl ^xl ^1' generalized 

coordinates^ ^ and S, and'' are thus the required forms of 
angular velocities and accelerations. 
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We proceed now to find displacement vector P. This is 
identical to X, U. 6 ). Replace A in ( 6 . 6 ) by W (Eq. 
(5.31)). We find 


(7.62) 

p /-lx, x„ 4 A\ f. AIa ^ Aa, + VV 4 R 

Expansion of (7.62) yields 

^ - \ - X 2 1^*63) 

^7 ^ ^ I lo 

p2 t A vv + s i o p ^ 


In the derivation of (7.63) from (7.62), we used (5.83) 
for V 9 emd (5.84) for W 9 to express (7.63) as a function 
of idv and^dw (See (5*75) and (5.76) for the expression 
for A^r auid 4 w) . 

To obtain (7.63) we neglected the elastic extension u* 
and third order products of small quantities involving Ai , 
4 ^ 2 , Oe, yiQ, and zj^q. Displacements v and w are second 
order products of elastic variables v, w, and Oe (see 
(5.75) and (5.76)). 

The retention of ^v and Av in (7.63) is subsequently 
justified by showing that they yield terms in the 
generalized force Qj comparable in magnitude to the 
terms contributed by the other terms in (7.63). 

Differentiation of (7.63) yields 
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(7.64) 

^ j A I 2. 

- V' t (Iw t C ^ ^ S -3 ^ 

^ 7 . - W i O^Vv' i St ''I 10 ^fj I o 


V 





I ^ 


(7.65) 


f ^ - V' ^ A V + ^ 9 I ^ ^ ^ 

f • * * 

T 2 ^ - w +- A w t' f V. - 2- j Q 


To complete the reduction of P^s# Pys# end P 25 end their 
derivetives to useful forms we require expressions for 
j\l, A 2 , end S end their time derivetives in terms of the 
coordinates w, v, and and the non-elastic pitch angles. 
Before doing this* we state the following definitions for 
pitch angles. 


-- &c + 

Q - G ^ G^(^ 


(7.66) 

(7.67) 

(7.68) 


where 





(7.69) 


Expressions for end A 2 derive from (5.60) end (5.63), 
respectively. 


114 

PAGE 


tpa ft 



Sikorsky Pircraft 


OOAMAATION 


REPORT NO. SER-S0912 


oiv«iON cm uNmo 



A, - - SqV ' * 


n.7o) 


The trlgonomentrlc functions required in (7.63) to (7.65) 
and (7.70) are 


^6 - ‘ 

'■ *■ 

where 0 ^ ^el considered small ezu>ugh to justify 
the expansions in (7.71). 


Substitution of (7.71) in (7.70), differentiation of 
(7.70), and neglect of products of quantities Involving 
V, w' , ^el assembly with the corresponding 

approximation to (7.70) yields 


/ 



- ^ 



v' 


) " 


, 1 

a ' 

C 9 ^ V 

. . 1 

^ * 

i 


# * 

A. -■ 

C-e.v'f 

i 

» e 

v^/ 


(7.72) 


Derivatives of trigoxrametrlc functions in (7.64) and (7.65) 

are 
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where 





- GO, 

, . ^ e- 





(7.73) 


(7.74) 


To obtain (7.73), we neglected products involving^, B, and 
^el* 

Equations (7.72) to (7.74)are the expressions which are 
substituted in the expressions for PxSf PySr and Px 5 and 
their derivatives, (7.63) to (7.65). The'^approxiiaatiorb 
(7.72) to (7.74) will be justified %d\en generalizec t • -c 
Qj is form^ by showing that higher order approximations 
would contribute negligible terms to Qj. 

At this point we have derived all tlie information needed 
for a \rorlcing form for the acceleration contribution, 
d^P/dt^, and with this done we are iu a position to 
present in the next section working forms for the to';al 
acceleration of the blade pi«rticle, defined by (7.4). 
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7.1 VTaions of Blad< Acc«lT«tlon Expression 

Pour different versions ere used for the blade particle 
acceleration expressions with a mind to easing and system* 
atizlng the subsequent derivations. The first version 
eases the evaluation of section integrals required in the 
expressions for the inertial loads Px§r PySr Pz§r qx$» qyS* 
and qz 5 * The second version provides an ^explicit notation 
defining the elastic displacment and pitch angle, 
orders of magnitudes of the terms, to aid the derivation of 
approximations. The third version separates terms into 
radially independent and radially dependent members to 
aid the evaluation of modal integrals in Q^. The fourth 
version facilitates the distinguishing of blade modal 
acceleration terms from modal displacm&ent emd velocity 
terms to assist the formation of the coefficient, s-tu, of 
the modal acceleration q']^ in (4.14) *' 

a) Section Property Version of Acceleration 

To facilitate the evaluation of section properties, we 
express the acceleration components in the X 5 syst«a in 
terms of coefficients multiplying the section coordinates 
yiO and *10 



(7.75) 


To derive the coefficients mb utilize (7.4) which 
reduces to 

X- o ^- ) p)xf- (7.76) 



following substitution of (7.54). Replacement of P^e, 
Py 5 , and Pz 5 # and their derivatives in (7.55) by (7.63) 
to (7.65), substitution of the resulting expressions for 


»• RPV D 
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(ao)x 5 » (aD>v5» (««> x5 (7.76), and oomparison of 

(7.75) and (7.76) yialas the following expressions for 
the coefficients in (7.75). 


(7.77) 




iJ ^ s 


t • * * 

0 X - ~^z ' ^ 0 ^ ^ ^ ':S‘y 

f ^ 3 \ 

^ f) X ■ \ - 

, -A, ^ 

-i- j I 

/\,-- ('=‘»,,-)cj ' S':) 

+ V i [\ v - 3 )^ -( W t A w ) f ^ ^ V ) 

' (^1-; 'x^- * ‘'-' l^.)(\- * 

<- X,' ^ '•“•^5- ‘~'j ,')(''^^^ 


^a rw n 



! 


I 
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a 


I ’ ^ 

b " 5^7 


- 

*- ) 

^<1. 

■/ ("-tox^ 

- -h u*., 

. 




-5^ - 3 ‘^ZyX, 

- a 












Sg ■/• (^- <0;j 

^ Co z . 
3- 



Ar 


r 





^w+av^;' f a -t- Av ) 


y* ^ 



)(w 


" 5q 

j V a K 5- 





- i- J. ^ Cjy 

r . 

^‘^^r3Aa 






‘"'sOcg - 

*■ 




0 » 

- ^9 f 

, - a cj ^ ^ s 

e 






coxj <-co^^ 
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We neglect the higher order displecenents v and w and 
their derivatives in (7.77) in all coefficients except 
Ay and Ax* This approxination is justified in the next 
chapter by showing that in the inertia loads P„R» Pxs# 

<lx3» Qy5f only Ay and Ax need to be approximates to 

this order tc ooi^>ly mth our approximations to the 
loading equilibrium or modal equations. 

b) Order of Magnitude Version for Acceleration 

To present a version of acceleration specifying the orders 
of magnitudes of msmbers in the acceleration expressions 
in terms of elastic displacements or pitch angle, 0 
define the coefficients in (7.75) as 


(7.78) 






I- rt 




0 , . 6 /') 

Cx : 

. 13 

A, = A/'Va, 

13 ^ , / 3 /-^ 
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Superscript (0) indicates that the coefficient is of 
zcuroeth order in tenas of elastic displacements or angle 
" that is, such a coefficient contains no elastic 
displacements w, v, 9m, or Superscript (1) indicates 

a coefficient with a first order dependence on w, v, 0^, 
and Superscript (2) indicates a term %d>ich contains 

products of elastic displacmnents w and v, and 0^, 

Si^titution in (7.77) of the ^proximations to s^, and 
C 0 their time derivatives defined in (7.71) and (7.73), 

and comparison of the resulting coefficients for Ax* 
etc wi^ the definition (7.78) yields the desired expressions 
for , Ax'^^ etc in (7.78). In these expressions Ai, 

A 2 rid their derivatives are to be approximated by (7.72), 

2 uid (7.69) defines ^e^ * obtain 


(7.79) 




A 


f5 


0 


(o) 

r 

cO 


- a V' i- ^ w' 


V 

t fco 

3 Jr *> / 


c 

( 

9c 

' • m 

Aj^ t 2 5£p^ 


^ CJ» 


3 f ^ 




\A/ 


’0, 


c 

c 


, To) 


X 

(>) 




- A , 31 co^ 5e^ 

a 
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The second order character of and A 2 ^^^ is due to 

the second order character of id v ^md which Involve 
products of w, V, and (see 5.75 and 5.76). 


c) Modal Integral Version of Acceleration 

To introduce a version of acceleration appropriate to 
the evaluation of modal integrals we express the generalized 
coordinates in (7.79) as modal sums, and then employ a 
convention distinguishing the radial dependences of the 
terms. 


The modal transformations substituted in (7.79) are 

/3--A 7c' 

S - S', 


(7.80) 


. 6 )^ (7.81) 

with the repeated suffix i indicating summation on the 
number of modes, M. 


Derivatives of these modal sums required in (7.79) are 


S, s 


A 7c 

■■ 7c- c <57 7 : 



-- 6^., A 6), 






?c' 




(7.82) 


(7.83) 
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Substitution of (7.81) and (7.83) in the expressions for 
ng , eg and their derivatives, (7.71) to (7.74), 

yields 


^ Cq 

\ wv '7J 

'"c- t 5 (J 

# * ^ 

A I ”* 1 / ' o ^ 

^ f: f <.|J^ V 7 


/\ 


V'. 'a 


<-■ ?t- o . 

’e ,e^. 

c g ^ B, f 

S^3 

H ^ '(^t ^ 7c ) 


(7.84) 


We also require the time derivatives of 4v and .Aw for 
substitution in and A, (2) in (7.79). From (5.75) 

and (5.76) ^ 
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(7. 85) 


(7.86) 


a 




A 


%■ k-fl -(/^Y^ j . 


(7.87) 


(7.88) 


Differentiation of (7.85) and 

A^> - ft ( j, cj^ , ) 

Av - Av-,H('f,?,, ( Sy; Y* 

A w- = A vv J (y^. ^ ) 

A - A lA,'y ( <^^. <J^4 -;i y . 9/^ / j 
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Thtt occurr«nott of two repeated suffixes, i and k, in 
(7.85), (7.87), and (7.89) indicates a double sunmation 
on the modes. 


We illustrate now the expansion of typical terms in 

(7.79) , to exemplify the derivation a version of the 
accelerations separating radially dependent terms from 
radially independent terms in the acceleration ooefficients 

(7.79) . 


Substitute as required in (7.79) the modal sums and their 
derivatives, (7.81) and (7.83), and Ai, ^ 2 , sg, c^ and 
their derivatives from (7.84), and introduce below a set 
of radially independent coefficients. We obtain, for 
example. 



(7.90) 




(7.91) 


V', % % -(^Xf ^ V <’•«> 






(7.93) 
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( jc ?« ' % 9 ,< 


+ (‘"‘'x,-* ‘■''2v‘''V) 7> 

^ Aw^ t Av^^ 


(7.94) 


(7.95) 




0 .) 

b 



^'c ) 

''c ' %, 


)Ce, 

* %,..% )?c 




>-‘^z/)( i9^^y 1 

^ ) 5e, 



' 


- < 6 ). ' ^// ' 9 ;. ♦- Q^, 




(7.97) 
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Comparison of (7.90) and (’’*91) yields the radially 
independent coefficients aol and Similar comparisons 

yields the other ooefficients. We obtain 


(7.98) 


) ■j y ^ 3 "I j' 

^ H.; ■►"/f )'?c 

-3 ? c f k ' ^ 5- ‘^'•V ) 

- -4 ’'^4 
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The symbolism for the coefficients aQl^ (7.99) 

indicates the origin of these coefficients, and may be 
used to trace the source of a term and to attach a 
physical intepretation >’:o it. A superscript x such as 
in aoi^ indicates a derivation frcxn the x component of 
acceleration, a^, in the X 5 system. Likewise, super- 
scripts y and z indicate ay and derivations. A base 
'a* indicates an A term origin in the accelerations, 
(7.75). A base 'b' would indicate a B term origin, - 
and a base 'c* - a c term origin. The first subscript 
in apj^^, etc in (7.100) indicates the order of magnitude 
origin of the tera. Thus aoi^ derives from A^t®^ , an* 
derives from Ax'^^ , and a 2 X^ derives from Ay(^). 

The occurrence of integer subscripts i or k in any term 
in (7.99) defines the summation treatment required. 
Generally, no i or k subscript indicates no summation. 
Thus the terms cqI^ s^, CQ 2 y cffr in the expression for 

in (7.99) axa not sums. The occurrence of a single 
i*^ subscript rewires a single summation. Thus the term 
an* W|^ in A]j^(i) in (7.99) is a single sum (Note that a 
axx^ contains a suffix i and thus the product a^x^ W£ 
complies with the repeated suffix summation convention) . 
Similarly, aox^^^Vj^k expression for Ay^^' is a 

double sum with i and k repeated in a2X^ (See 7.100) . 


d) Version of Acyleration ^pression Separating Modal 
Acceleration ^ from g and q Terms 

In this version of the acceleration expressions we 
separate the blade modal acceleration rerms involving q, 
from the blade modal displacement and velocity terms 
involving q and q in acceleration coefficients aoi* etc 
in (7.100). This version facilitates the formation of 
the coefficient s-;)^ of the blade modal acceleration, 
q , in (4.14). 


The only sources of blade modal acceleration terms 
coefficients aox* etc, (7.100), are^^xS' ^5 **^z5» <Ii 
and 0^. Enqploy (7.60) for 4 »/.x 5 ,. 4 i^y 5 , auid «X>z 5 , substitute 
the modal transformations 4 , .,S , (7.82) in (7.60), 

introduce the expression for 0^, and underline the modal 
acceleration terms. We obtain 


6J 


^0' 








(7.101) 
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' p’k'IK 

^3 )k 


- ^Cl, '-^K 


f^l S Cf ^/S ^ b ^ ^ 

^ (^'H)g fa^ S\ ) 9 ^ 

The symbol is the Kronecker delta and should not 
be confused with the lead. angle component of the mode shape , 
The expression for derives from the expression for 
St ostablishcid subsequently in Chapter (12). 

Remove the modal acceleration terms from (7.101). We 
obtain the following residues, denoted with asterisks. 


cO , 

^ / 3 c^ir' 


(7.102) 


* * jk 

/ L - o 


. . 0 ) : . . 


^ fc " f^tS ^ ^ ^ 


/s f Oy^ 


Equation (7.101) may be expressed as 
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(7.103) 
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To develop an explicit notation for the^ residues of the 
acceleration coefficients deprived of q* terms/ and for 
the removed terms themselves which are q’ dependent/ we 
re*define the acceleration coefficients (7.100)/ in terms 
of residues £md q - dependent terms. 


We express the coefficients in (7.100) as 


(X 
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a. - of ^ )C 


(7.100 
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(7.105) 


All the coefficients in (7.100), except the four listed 
In (7.105) Eure lilce (7.104), and may be generalized to 


CL 




(7.106) 


mn - 

{ X X V • 

^tnn - 

r ^ X ^ ^ • 


In (7.106) X, y, z in (7.100) substitute for X to yield 
from (7.106) any required element in (7.104). The Integer 
subscripts m and n are taken from the subscripts in (7.100) 
to generate the elements in (7.104), with exceptions not 
complying with this rule listed in (7.105) . 


We illustrate now the derivation of typical terms <afmp 
>^mn*» Vmn* (7.104) and we al 

iXlustr^e the derivation of terms in (7.105). 


X 

# 

SO 


We substitute in the coefficients AqI^ (7.100), the 

terms iA'xs, ^v5 from (7.103), as required. Typically, 
obtain fr^ (7.100) 


we 


’^‘^ 5 ,- i-(‘^ l( 'Pk fk*- 


(7.107) 


r cX 


1/ 




(7.108) 
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Comparison of tarms in (7.107) and (7.108) yields 


(7.109) 




Similarly, vfe find 




= A 5 - 7 


yielding 


1 c 


p/3 ^ 6 


CJy^' 4- ^^'z\)7c 


5 


The second coefficient, ^ 13 ^ » comes from 


/ 


7 . 


' K Ik 


yielding 


s 


Lk 


(7.110) 


(7.111) 

(7.112) 


(7.113) 


(7.114) 
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Detemlnatlon of tho ooefficientt in (7.105) is 
exemplified by 


a 


21 


I- ‘in lilK 

L*ik pcp S' <• 5/3 — 1 


(7.115) 

(7.116) 


Comparisons of terms in (7.115) and (7.116) yields 

(7.117) 
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(7.118) 




^ 


The notation k <4> 1 in (7.118) means that k must be 
replaced by the new dummy suffix 1 in and in 
multiplication of qjc'^'l with terms like dwijc, subsequently 
performed, reminds us to change the k in^w^v to /. The 
reason for the replacement k <~1 in (7.118) is to prevent 
the occurrence of a triple subscript k in (7.118), which 
has no meaning in our repeated suffix convention. 

Equations (7.119) and (7.120) given below list all 
coefficients o^mn^» '^mn*r dmn^» >^mn*» ^mn* 

derived from the acceleration confidents in the above 
manner. 
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We recall that and in (7.119) are the 

residues of «mn » bmn^ (7.100), after removal 
of the blade modal acceleration terms q. Consequently, 
it is reasonable to expect that the substitution 


a 

b 


»vi A 
X 

/n n 



C 


X 

/II /t 






(7.121) 


in Qj^, formed in the next section, will yield , which 
is the residue of Qj“ after removal of the q terms. 


Similarly, the substitution 


a 


* 


X 


c 


X 

tn K 


yn r- 


(7.122) 


in the expression for Qj° should yield Sjj^. 

The validities of substitutions (7.121) and (7.122) are 
demonstrated in the next chapter, and these substitutions 
are thereafter employed to facilitate the derivations 
of Qj^ emd sjk* 
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7.2 Ordyre of Magnitudes of Acceleration Coafficiants for 
Small Hub Accelerations 

V 

In this section, we illustrate, for small hub accelerations, 
the procedure used to determine the magnitudes of the 
acceleration coefficients in (7.119}< and (7,120). This is 
followed in Section 8.4 with a description of the procedure 
for neglecting small terms contributing to for small 
hub accelerations. The elimination of these*^ terms effects a 
substantial reduction in the number of significant terms 
retained in considerably simplifying the algebra. 

As a preliminary %re state that € will designate any small 
quantity, while unity designates a non-small quantity. 

To nudce sure that the symbol s , which uniformly designates 
the orders of the different parameters, assumes comparable 
values 2 UK>ng the parameters, we non-dlmensionalize the 
modal equations. In the following section, we assume that 
all quantitites have been non-dlmensionalized, according 
to the method of Chapter (12) , without introducing a new 
explicit notation for non-dimensional quantities. 


The assumption of small hub accelerations requires the 
following axis quar^-ities to be small. 

•> z c (e ) (7.123) 

* » 

In addition, it is convenient to recall here that the 
following are also small. 


, S ,w, , \/, 

(Pf , vJ ^ V' 

\</\ ) ^ J t 


~ 0("6) 

r O ^ ^ I 


(7.124) 

(7.125) 

(7.126) 


w ft 
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W« note that since wi » w and - v + Sr, and ^ 

and S are of 0(€), then and are of 0 (e ) in (7.124) 

Dlaplaconents 0^, w, and v in (7.125) are expressed as 
elastic dlsplacMient orders in contrast to an ^ -order to 
achieve the specif icnese required by our assumptions 
distinguishing the orders of magnitudes of 0^, w, ^md v 
from f . These assumptions are invoked in the next 
chapter to derive Qj“. 

We illustrate now the derivation of the orders of 
magnitudes of typical acceleration coefficients in 
(7.119) and (7.123). 

From (7.119) 


^ (7.127) 

(7.128) 


Employ (7.27) to obtain (aos)x5' (7*39) to obtain gx5* and 
(7.59) to obtain *Jy 5 ando>25. find 









(7.129) 

(7.130) 




.n 


if 


cJ 


f a 






-h 




^-■2y/3 


CO 


z 


f 


a 

3 / 



f ci. 






(7.131) 

(7.132) 


Equation (7.19) with (7.24) yields (aos)xl' ^^QS^yl 
(ao5)zl* Invoking the order assumptions (7.123) and 
(7.1261, we find 
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(7.133) 




Similarly from (7.37), (7.39), and assumption (7.126), 
we find 


J K, 






0 (c-) 

0 (0 

0 C^-) 


(7.134) 


Employ assumption (7.124) and the expressions for the 
direction cosines (7.29) to obtain the orders of 
magnitades of these direction cosines. We find 


(X 




I 

o ( 1 

e 


1 

I 


e 

O ) 

I 


(7.135) 


find 


(7.133) to 

(7.135) 

in (7.129) to (7.132) 


- c 



- 0 

CO 


r C 

CO 


r O 

(0 


We 


(7.136) 
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Substitution of (7.136) in (7.127) and (7.128) yields 

(7.137) 


^0/ " 


o ^ 


(0 

0 Cl ) 


(7.138) 


Similarly, we find 


oL 


II 


- C 


; (« c ) 


(7.139) 


(7.140) 


and so on. 


From (7.120) we find 



(7.141) 


(7.142) 
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The orders of all terns in (7.119) and (7.120) are 
similarly derived. 

It is necessary to )ceep the modal displacement qi in the 
order statement to help us distinguish between terms of 
0(£) and terms of 0(^q), 0(w), emd 0(v) when Qj^ is 
formed. For example, when a term O(q^) oomblnes with a 
term involving the resulting order of the product is 
0Q because « ^i<3ir while a product ^1 " S 

is of order e . Because terms of 0(fi*^), o(w), and 0(v) 
are approximated differently from terms of 0(6)* we have 
to retain q^ in the order expressions until a point where 
products li)ce ^iqi* Wj^qj^, 0141* otc are formed in the 
expression. 

To emphasize the physical Implications of our order 
statements we list a few examples from those given above, 
and revert to the interpretation that all quantities 
presented below eure dimensional. Thus we have 


(7.143) 



- o(e') 


= 0 (^) 

/Jl"" 

-- 0 (t ) 

<!’ jlK 

- 0 (e) 

^ jjl 

= 0 

u; j ^ /ji 

= cO) 


-- 0 ( c-) 

ox /jl 

^ 0 ro 
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8 , Inertial Generalized Force Q jP 

We derive In this chapter a working form for the Inertia 
contribution to generalize. 4 force, Qj^. 

We begin by deriving expressions for Inertial loads 
Py?'_Pz5' ^x 5 ' terms of accelerations in 

tne X 5 system, il^.75). Tne resulting expressions Involve 
the coefficients Ay, Az etc In (7.75) and section 
properties like c.g. offset and radii o.^ gyrations. The 
inertial loads axe then substituted In the ganerallzc ' 
force expressions , Fy, and F 3 and Integrated to form 
Qj“^, and Qj®^, (4735). The different acceleration 
inventions of tne previous chapter participate to effect 
approximations, develop modal Integrals, and to separate 
dependent and q Independent terms to form Qj^ and sjk* 

8.1 Expressions for Inertia Loads 

Here we derive expressions for the Inertia loads per unit 
length of blade, p^J; PyS» Pz5f ‘Jy?' terma 

of acceleration coefficients Ay A 2 etc of the system 
and section properties, comprising c.g. offset and^sectlon 
radii of gyrations. 

Figure 17 Illustrates the positive conventions for the 
Inertia loads resolved to the Xxo_dlrections. In terms 
of accelerations resolved to the x^g directions, these 
Xxo Inertia loads are 


PX/o 

(- ,^pd ^ 






(S.2) 

c— 

- 




1 I/O ■ 
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Figure 17. Positive Conventions for Inertia Loads 
Employed in Multi-Blade Rotor Analysis. 
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To derive X^ system loads in terms of X 5 system accelera- 
tions, we employ the following rotation transformations 
relating free vectors X 5 and x^o* 


(8.3) 


(8.4) 


Free vector X 5 _is selected from among P 5 , q 5 ^ emd a 5 , and 
corresponding X^o vectors are p.., q^^g, and aj^g, with 
components 




(8.5) 


t> > f 3 ^ % y- 


a 




a 


2 . 







Application of (8.3) to P 5 and pig, substitution of (8.1) 
fojf_PlO substitution of (8.4) to express a^o in terms 

of a 5 yields the following set of results. 


r -- 

- (- Agp a y J A 

X' P d A 

( 8 . 6 ) 


160 

PAGE 



1 


r 


f 


( 


Sikorsky Pircraft 


u 

iumcm. 

p 


REPORT NO. SER-50912 




The expanded form of (8.6) is 

«.x- 


5 

q//-) . f o. M ) 

r. 


(8.7) 


t? j 

To obtain qs in terms of a^, we first write (8.2) in the 
form 


7,0 = Sp 


% ( 


o 

-t. 


^/C 
/O o 

^ I O 


o 

o 


) 


( 8 . 8 ) 


/<- 


Following a similar procedure to that used to derive (8.7), 
we obtain 


7y ^<^7/0 

‘J,o 0 


0 !3 /• 0 




(8.9) 


This reduces to 

?X^ ■ ,r_ ^ ^ J ~7 a q 


-/4 


( 8 . 10 ) 




( 8 . 11 ) 
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-T 






(8.X2) 


We neglected third order products of elastic variables and 
fourth order products of small quantities to derive 
(8.10), and we neglected second order products of elastic 
variables and third order products of small quantities to 
obtain qy^, (8.11), and q^e, (8.12). These approximations 
are discussed in more detail later in this chapter. 

To reduce (8.7) and (8.10) to (8.12) to forms involving the 
acceleration coefficients Ax, Ay, etc, we substitute (7.75). 
Assume inertial symmetry of each section and introduce the 
following definitions. 




(8.13) 


ir 


5--=^ X, 




We obtain 

O 

Fx.' 


- ty> 




(8. In) 


o 

S' 






(8.15) 
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^ ‘3c' - [(^ Pf 















(8. 

.16) 

— 

» 







f ( I3^0g + 65 





C <r - 



(8, 

17) 

Shears 

derive from (6.71) 






- AT ^ . 

,-W f 
) ' 

(8. 

,18) 

- X 


- - ^ tn 

c- 

vJ 



'*?'X 

V /) 

" ' 4 - r^ 2 . ^ 





Loads and tensions in the Xiq system derive from (6.60) and 
(6.61) . 


/ y 


^ 'P^'? 


(8.19) 

F.; ^ ^ ^ 

k ‘ ('2/ 

(8.20) 


Substitute in (8.14) to (8.17) the expansions for Ay, 
A 2 etc defined by (7.78), and employ ^ 
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- + Ge, 5 ^ 

- ( 9 |_ r- ( 9 ,^ 


( 8 . 21 ) 


We obtain the approximations 


?X,- fss^rc<5^ -J.u'hAf'l'- 


( 8 . 22 ) 


+ nS°(<=s, - 5s., )fC ‘'Ve^ . 6c, Ce, )f CV'4,)fc 






(8.23) 




^kC 9 , k? 

'® ’- J fo 


(8.24) 


■. ‘ - H"«"(-.,e-, >»,) . s,% . «;''x,iv 


^ R.^°) 


^ ^ X S 0 ^ 
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with ^ 2 ^ and A 2 approximated by (7.72) 

A,-- v' +C6, w' 

A;l"- *^64 V ' <• w ■ 

Corresponding approximations to P„ 5 » F^, (8.18), and 

PxlO' ^xl 8 » ^®* 20 ) are ^ 


f ^ 

-- D 

■ 




SM 

■r-[ 


Afi A 




(o) 


^•y frx (o) 


/o. 


(') 


(8.27) 





and 




D 

10 


- - K>T 




f V t, 


(8.28) 


r; 


o 

/O' - 


r\[ 

r L 


J 


, ^ 


/ ^ ) 


/I /] <^'7 f ,s '"' r> ^ 'y ^ 


; //7 • ,7 <'V 

^ ^ ^ ^ 2 . 




C.W 


■ 


(8.29) 
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In (8.22) to (8.29) our approximation to all loadings 
except are accurate to first order elastic terms w, v, 
and 0 ^, and second order products of small quantities 
involving w, v, 0 ^, y^ocg' ^^ylO ^ *^zlO* Loading q^^§ 
is approximated to a higher order than the other loadings 
and is accurate to second order products of elastic terms 
and third order products of small qua. :ities. We show in 
Section 8.4 ti at these approximations when substituted in 
the modal equations are sufficient to satisfy our 
approximations to the modal equations. 


8.2 Derivation of Typical Terms in Q j 

In this section we illustrate the derivation of typical 
terms contributing to the generalized force, Qj, to 
exemplify the procedure used to derive a working form for 
Qj, and to introduce the modal integral notation. 


We base the formation of Qj on (4.36) and (4.35) 


(8.30) 






r' O ;l 

V\// ' 

J 


r* . ,-r o . 


(8.31) 


Substitute in (8.31) the expressions (6.63) to (6.65) for 
to P 3 ^. We obtain 


‘ D 




-((r- ) 
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- ‘-■fl V " 


t V " 5,, e 


I-' “ 

t/ ^ QW 


C’' ^ 

J 




/ 0 >Ol 


d'if 


"f w- 


[\e. 


0 y 


p C» 




J)/. 


' C V/ 

-(El.. EI^M.Vf„,iJ.-;„-j!:„;jj. 


J. w. f,»; ,(vT,,“)'. p.» . r-fw/) 


(8.32) 
/s / \ / 


■f- tyi vv'. 


Cvv' ^ 


V ( '^A ^ Ac f ^ 


C w' 


(8.33) 


J -'j [-1^5 V(v ey 3- )* ..;)-»- JlV,+).v 

- l/'Q *^A '"x ° *■ A Vv^w ^ 0 , ^ ) 


C\A/ 


' &e -e (- ^',3 ''x 3 ° t A ep ^ y .v/^ 


® <• Vi/ 


(8.34) 
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We illustJrate now the expansions of typical terms contri- 
buting to to Qj“3, (8.32) to (8.34). Use (8.23) for 

qx5° in (8.32T» employ B’* 0^ ^t ^ expand the pitch 
angle, and neglect third order products of elastic variables, 
w, V, and and fourth order products of small quantities. 
We find 

c 


-i >. V -j ^ ^ 6., Ce, V 

-[Az (ce, • '9-., - J 9efCe,*y,^')^ 

^ e-, ♦ 8 , '-'> 5 ,^ ■ ) _ 6 /'^Ce, 

^ ^ ^ ■ 'A.-, 5t.. ) + C>|‘Ve^ ^ f '■^f, ) »- Cz ^ 




^oc "* ^c>> k 5^^ j 


(8.35) 


In (8.35) and A 2 ere approximated by (7.72) (also 
duplicated as 8.26). 

Employing (7.99) .to expand the acceleration coefficients 
Ay , Ay'*^, Ay' , etc. in terms of radially dependent 
and radially independent terms, and employing 


^ -- ^ 

- ^ 7 c' 


(8.36) 


where required, we find from (8.35) the typical development 
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^ ^ &<p, - i 


h T 






(8.37) 


Denote modal integrals originating from Qj^l by 
etc. We find that in (8.37) yields 

- 4 ^ Tt ?/< J 5s^ ^V 

o/'4 q (‘) , <). ^ 




Q,% ^9^ q(')^ 

^ <9^ cp/'J, (j 0) 


gf = r^'^irK^y ^<9c cCr 


’ c 

(8.38) 

(8.39) 



r- 

of -- 



0 

, (0 
9(, •- 



c^>- 

'5 
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The contribution to (8.37) from is 

0 


1 

<? 


- CL 


oa 


c; 


V 


03l 


Q 


Ci) ^ 

e ^ o a 


q (p^‘^ 

7c 


(0 /'. ^ c> 




40) 



(8.41) 


Similarly, we find 

5 'i&j In A .j'' V^e c '9'f, )"5,0(. , <ii' - 

C/ ^ ^ 

S^ ^j ^ ^ ^'c ) (^0 c ^ ^ ^ ^ 

0 u 

(8.42) 
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Th« contribution to (8*42) from is 




Q 

Q 


n ^ lb 

0 ] 

lb 


r ^^(90 


o 

r' 






(8.43) 


(3.44) 


I 


The contribution to (8.42) from ax2^ is 

e ^a.'’ 

J 


7i 




* -,t o/;' 






(8.45) 


0 ^''’ -- 


/t 

J Oj r»w, 

'' -A 


(^*^r 

'■ 

J, e --"^c ^Cc b.-o./'- 

0 

^/8 ^ , 

1 yz-c-s 


Another characteristic term is 




ci) r * 
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/O ^0 

^7 


/ 




t) V/ 

-- J '‘0 6 >j,0^ 


(8.47) 

(8.48) 

cwr 

J 


We,«valuate now the contribution to Q from - ( (T-P^s ) 
0L in (8.32), to illustrate a typical application 
of integration by parts to collectively differentiated 
terms. The contribution to Qj^i from this term is 

'/ ^ fT - " 

c 


( Qj (t- fx° ) Kfl e^‘ f Yt . F^°) kI ^4 


(3.49) 


The condition 

^A. 

T ~ ^ - C ^ Y- z (Ty. (8.50) 

was used to obtain (8.49). Replace in (8.49) Px5^ (3.2/), 
substitute ^e' * remove neglic -*^le terms. We 

find 
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^ f \ ^c' % 1 4 [s Y. . /; . a 4 ]j^ 

Inserting (7.99) for (8.51), we obtain 

(ji=e 


/a t 


o J ^ t / 

-(Oj e/ v/?j, ^ 

/?, cj i£ 




(8,51) 


(8.52) 

(8.53) 



parts. 
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BxamplM of tttrms contributing to Qj^ and (8.33) 

and (8.34), are considered below to illustrate the use of 
superscripts to denote the origin of modal integrals, and 
to prepare the ground for the examples in the next section. 


The contribution to from P-c® in (8.331 follows from 

a similar treatment to that used for the contributions. 




0 

VJ, iry, t- f ^ / 1 .-. ^ s 


<- [_bj Sg,^ f *■ 


^St hi ^9^ *^,X ■»■ k,^ &i ^9, 




Contributions from ag2* *nd bi 4 *, for example, are 

r»~r 

c/ 1' - n ^ C^/ 


■I ' ^/y '- ^/y 

jCu) _ 

~ J ~ W, irVir- <ih 

U) p 

" J - W, 

c J ^ 5)” 

The contribution to from mvj^ in (8.34) is 

r 


r' cir z cx.^ Q 

6' .3 ^ 


(8.55) 


(8.56) 


(8.57) 
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(8.58) 


The exanqples (8.55) to (8.58) show the use of superscripts 
2 and 3 to denote integrals originating from and 

To aid the identification of the physical origins of terms, 
subscripts attached to modal integrals increase progressively 
from a value of one upwards as one progresses among the 
terms composing in (8.32). The first terms in the 

integrand of (8.32) contribute the lowest subscripts in 
the modal integrals, terms succeeding these contribute 
higher subscripts, and the final terms yield the highest 
subscripts. Equations (8.35) to (8.48) illustrate this 
convention. To preserve this sequence, modal integrals 
reducing to the same expressions are not given the same 
sxibacript if they originate from differcmt places in (8.32). 
Exactly the same .conventions apply to modal integrals 
derived from and (8.33) and (8.34), 

Table (1) relates the physical sources of modal integrals 
in the integrands of (8.32) to (8.34) to the subscripts 
attached to the modal integrals, and may be used to identify 
the sources of terms contributing to Qj^. 

As always in this report, the interpretation of terms 
involving repeated suffixes in (8.38), (8.40), (8.43), 

(8.45), and (8.47) must comply with the repeated suffix 
convention for summations. It is necessary to Iteep in 
mi^ that suffixes may be burled in terms ajx^f ^ 21 ^ • 
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//V nr (ea .?.?J 


So^^(/?/pr\ 
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(Vi 
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- <^,J I (- 
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4,V 
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33-36 

37 

So 

S/ -/34 
/ 3S " /So 
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8.3 Inertia Force Components Q j ^* and S jj^ 

In this section, we separate the terms composing Qj^ into 
a component Qj idependent of modal acceleration, and 
a component s^]c%# isolating the modal acceleration, and 
we list here all such terms in and Sj]^. The derivation 

is specialized for small hub accelerations emd this assump- 
tion effects a substantial reduction in the number of 
significant terms retained in our equations. We present 
also a classification dividing components of and sjK 

into terms not involving modal sums, terms involving single 
modal sums, and terms involving double modal sums, to 
facilitate the programming of Qj^* and Sj]^. 


Our aim is to form the elements (4.42) and (4.15) in the 
modal equation, which are 




D * 


s- 

"J K 

with 


- ( aJ (3J 

^ CL • 4 s + S 

J f< JK ^ J K 


(8.59) 


(8.60) 


' *0 ^ ^ ^ 3 /o ) ©j Wi^ + v,^ Vi^j a', 


ar 

(8.61) 


The generalized mass a^K is known from the solution for 
the normal modes, and the other elements in (8.59) imd 
(8.60) are derived here. 

We proceed now to illustrate the separation of Qj^* and Sj]^, 
drawing on the examples of Section (8.2) for the ''derivation 
of typical contributions to QjO* and Sjv. The basic step 
is the substitution of acceleration coefficients (7.104) 
to (7.106) in the ^nas composing Qj^, to effect this 
separation into Qj° and Sj]^. 
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□•note with an arrow a contribution to QjO or from a 
subset of the elements making up Qj^ or Sj]^. We may write 
the contribution to from the terms in (8.40) 


as 


Oj - 




V O. 


0^ w., i G^Q,' t. ^ 
0 ) 


(0 

'3 


0.? t ^ 7 


^,o' t- -o'! Gt C?^'^ 

)A 






(8.62) 


Employ (7.104) or (7.106) to separate the acceleration 
coefficient ap2^ into a term independent of and a term 
dependent on q]^, that is, 

. (8.83) 


Insert (8.63) in (8.62). We obtain 


oa 

-I- rX 

oa 

— y) 


Q 


0 ) 




(■) 


f o< 


(0 


(•> ■■ 0) 


9\ Q ‘"-/o ^ Q Q n9i ■ r^‘ '> 

''?a ?.r -!%x ?K 7^ 9i<^ 


(8.64) 

The notation k « 1 is designed to prevent the occurrences 
of three k subscripts in (8.64), to preserve the summation 
convention, as explained in the remarks following (7.118). 

Comparison of (8.64) and (4.30), %diich is 


a 


Q. ‘ - C 


0 


) / 


'J 


- s 


o)- 

Jf< 7 


(8.65) 
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yields 




o ^ 


V- 


cA 


5 


J i< 


CA 


0 ) 

7 




o a 




<1- 

oa. I L cj 




OJ 


' o 


^ 7 


(0 


( 8 . 66 ) 


'o a 


6. CJ)f f 


^ p.l u 0 


('] 


(8.67) 

Equations (8.66) and (8.67) are illustrations confining 
the use of substitutions (7.121) and (7.122) to obtain 
QjD* and sjjc from Qj^, We oh 
tne contribution to Qi°^» (8. 


observe toat ao2^ 


Similarly, ao2^ 


( 8 . 66 ) . 
(8.62) , 


.62), will yield Qj , 

Q 2 ^ in the con^ibution to 

will yield the contribution to Sjj^tl), (8.67). In the 
subsequent illustrations, it is to be understood that 
substitutions (7.121) and (7.122) have been applied, 
without further preliminaries. 

We underline in (8.66) and (8.57), and in the subsequent 
derivations of this section, negligible terms for small hub 
accelerations to remind us that such terms make no 
contributions to our set of terms composing , presented 
below. Section (8.4) illustrates the procedure for 
establishing the magnitudes of terms contributing to Qj^ 
and sj]( for small hub accelerations. 

Employing similar procedures, we find the contributions 
from terms in (8.45) to be 


0 ^ 
of' C- 


b,. 0) 

/t 






(A 




n 


/a. 


?A 




( 8 . 68 ) 


5 




O), 


(8.69) 
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A «lmpli£ioatlon of (8.69) is achieved from the substitu- 
tion ^3^ « - ^{ 1 . (Eq. 7.120). Rewrite the nodal integrals 
in (8.46) as 




"■ %c cCa. 


(8.70) 


Substitute in (8.69) the result • -Sn^- We obtain 


' ^ C H "iJ '^/(l 




(8.71) 


(8.72) 


(8.73) 
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Th« propartics of th« Kronacker dalta raduoa tha singla 
and double suma in (8.71) to singla tarns and a singla sum, 
respectively, in (8.73). Tha progran exploits these 
properties of the Xronecker delta, wherever it occurs, by 
using forms like (8.73) instead of (8.71). Both forms 
eure listed in the sjk tables. 

A third example is the contribution from (8.47), which 
requires a slightly different treatment because the 
expansions of a21^ and S22^f (7.105), are exceptions to 
the usual expansions (7.104) or (7.106). Substitute (7.105) 
in (8.47). He find 

(8.74) 



0 ) 




With 


'-■'0 ■' - T ^ ^ '■ h , 

0 ) 


(8.76) 


Oir 


^JlO ^ j ) ^ J / ' 6^ . ^ /I v;' S 

we obtain from (8.75) 
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Aa indicatad above, the undcur lined terms are neglected 
for small hub accelerations. 

Proceeding now to the remaining exeunples considered in 
Section 8.2, we obtain the following contributions to 
Qj® and sjK* 

a) Contribution from (8.52) 



X' 







/ jc (j 


(8.78) 

(8.79) 


since 




- 0 . 


b) Contribution from (8.55) 


(5 

<^4 

^5S 

(8.80) 

<r C^) 


Q 

(8.81) 




(8.82) 

^JK <r~ 




- 

' Sck 

.. ) 
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c) Contribution from (8.57) 



(8.84) 


(8.85) 


These examples complete our illustrations of the 
formation of elements composing Qj^ and Sj]^. 

We present now the division of components of and Sj]^ 

into terms not involving modal sums, terms involving single 
modal sumSf and terms involving double modal sums, to 
facilitate the programming of and sjk* At the same 

time, we introduce a notation indioa ting‘d explicitly the 
subscript dependences of the terms composing Qj^ and Sj)( 
to indicate clearly to the programmer where subscripted*^ 
variables are required. 


Recall 

express 


that^ repeated suffixes indicate summations. 
Qj“ and Sj)( as the sums 


r 


U u 


^ <? 


0 ^- 


J 


f 






We 


( 8 . 86 ) 






(8.87) 
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( 8 . 88 ) 


(8.89) 

(8.90) 


(8.91) 
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(sj 



8 


(i) 


C 


J 


Table (2) Hats all the elements derived in the manner 
C'f the examples given aOsove, and tabulated according 
to definitions (8.88) and (8.91). 

Referring to these tables, we see that index m ranges 
over the number of elements in a column of Table (2) 
2 uid since m is repeated always, %re must sum on m in a 
column to evaluate any term in (8.88) and (8.91). The 
other indexes i, j, h range on tiie number of modes, M. 
Appendix 14.3 defines the modal Integrals required in 
Table (2) . 

We select now a few examples from these tables to 
illustrate the correspondences of elements 
» ®kmj » etc., in (8.88) 6und (8.91) 
with the elements derived in the illustrations, given 
above. 


a) Terms not involving modal sums 


n 

h 


0 ) 




- <?( 


y 


Q 


o) 


/<; ^ oj 0) \ 

I 73, J qO) 

b) Terms involving single rooual sums 


0l 


(8.92) 


(8.93) 


/ 

~ L ^ 

^ W 


(8.94) 

/3 

' A/c 

ij 



(8.95) 
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b) Terms involving double modal sums 



(8.96) 

(8.97) 



(8.98) 

[ 


Because (Sj]c^^M2 small hub accelerations is expressible 
as the more convenient single modal sum, (8.98), we do not 
require the form (8.97), emd its components are not listed 
in Table (2) . 
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table 2 continued 
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1) Verify the adequacy of our approximations to section 

loads Px5“» Py5°» <3x5°' ‘lyS®' cited earlier 

in tiis chapter. 

2) Illustrate the procedure for neglecting small terms 
contributing to Qj^ and Sjic for small hub accelerations. 

We osnsider first the approximation^ (8.22) to (8.29), to 
the inertia loads. We recall from Section 4.2 that the 
torsion modal equation (4.61) excludes third order products 
of elastic variables and fourth order products of small 
quantities. The flatwise- -Higewise modal equation, (4.65), 
excludes second order products of elastic variables, £md 
third order products of small quantities. We verify the 
approximations to Px5^# Py5^' etc. by inspection of (8.32) 
to (8.34). We see that qics^ in (8.32) must be approximated 
as (3.23) which includes second order products of elastic 
varial)les and third order products of small quantities. 

The occurrences of Vqys'^ end w'qzsO in (8.32) indicate 
that qy 5 ^ and qzs^ neeSL to include only first order elastic 
terms and second order products of small quantities, such 
as in (8.24) and (8.25). Similar inspections of the other 
terms in (8.32) to (8.34) justify the adequacies of all the 
approximations, (8.22) to (8.29), for Px5°» Py5*^» 

qx5°/ qy5°r qz5°» fx5°» Fy5°» ^25°' PxlO®* ^xlQ®* 
Inspection of (8.22) to (B.29) verifies also the adequacy 
of approximations (7.72) for and A 2 * Finally, approxi- 
mation (8.23) for qx5^ justifies the assertion of Section 
7.1(b) that only Ay and Az require expansions to second 
ordtfs (A^^) and Al^M in (7.78). All other coefficients 
^x '' ^x^' ' etc., in (7.78) occur as zeroeth and first 
order terms in_Px5^' ^5^' etc., (8.22) to (8.29). We 
recall that A^^' and a| 2) are functions of Av and^w 
(see 7.79), which are second order translations accompanying 
the rotation (see Section 5.3). Consequently, the 
presence of these terms in qx 5 ^» (8.23), shows that these 
second order displacements are comparable to other terms 
contributing to qx5D* This justifies the inclusion of^v 
and in our calculation of blade accelerations. 


We illustrate now ^e procedure for neglecting small terms 
contributing to and sj]( for bmII hub accelerations. 

We consider firsf^o^02^ ^i ^11^ (8.66). From Sectli 


(7.2) wj 


er rirsr 0f02-^ ^11 10. do;. rrwm aecucr 

Obtain the magnitudes of the components in ^ 
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07- 


- 

f- 

( L ^ n - 


o(e) 

o(t) 


(8.99) 


( 8 . 100 ) 




o 

' 3 

^ 0 ^ ) 

Combining (8.99) to (8.101), we find 

?c ^ 0 ^Oe) 


( 8 . 101 ) 


( 8 . 102 ) 


which according to the assumptions listed in Chapter 2 is 
negligible. Similarly, we find the magnitudes of the other 
terms in (8.66), and we deduce that all terms must be 
retained except those involving Qii'^' and 

To approximate the terms raa)cing up sj]c, we employ only the 
assumptions that blade flap and lead angles, ^ and S are 
small. We may restrict our approximations to these angles 
because these are the only motion parameters present in Sj^. 
By following a similar procedure to that used to establish 
the magnitudes of (8.102), we deduce from the product 
Sj]c^^)q]c that all terms in (8.67) should be retained 
except those involving and 0^2 

When hub induced accelerations, blade offset, and rigid 
blade emgular deflections are among the small ^anti ties 
considered, we approxj^ted the terns in (Qj°^ )q to O(c^), 
and the terms in (Qj°^*)o *hd (Qj^^*)o each ^0 ( 6‘^). 

Ly, we approximated the terms In (sj)c^^^)i 


similarly, w«* aipp£WAXiHat>«u «.u« 

0 {€*), and the terms in (sj)^^^) and ei)c(3)Toqk each 
to 0 ( € ^) • These approximations may be inferred to allo' 


kk to 


allow 
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somewhat larger hub and blade motions for systems modeled 
with rigid body flap and lead motions without pitch freedom. 
These orders, however, are not consistent with our 
approximations excluding from the general blade model other 
terms of oon^arable order, like those from cyclic pitch 
inputs. For this general model, the policy of retaining 
such high orders in the %iorking equations is sound provided 
applications are restricted to problems for which the high 
order terms are insignificant. It is well to note also 
that because these additional terms do not require 
summation on the modes, no significant penalty in program 
execution time is incurred by them. 
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9. Aerodynamic Generalized Force Q j^ 

The calculation of the aerodynamic excitation, Qj^ is 
organized to yield Qj^ from the calculated motions and 
tables of non-dimensional Ci, Cd/ and Cjoq/i^ characteristics 
as functions of section angle of attack/and section Mach 
number. The table look-up procedure is a consequence of 
blade element theory, which is assumed to apply to our 
model. 

The first section of this chapter contains a derivation of 
as an integral of section lifts, drags, and pitching 
m^ents. The second section of this chapter con^letes the 
information needed to evaluate numerically the integral 

We derive there expressions for flow velocity, section 
angle of attack, and Mach number which then yield the 
sections lifts, drags, and pitching moments in Qj^ from 
tables of these characteristics. 


Unsteady aerodynamic states eu:e represented by tables of 
section aerodynamic characteristics as functions of angle 
of attack 0 ( , and the puameters A and B of Ref. (8), 
proportional to ^ amd w , respectively. The time deriva- 
tives «< and M are formed by numerical differentiation in 
the program, and are not discussed further in this chapter. 


In addition to the restrictions imposed by blade element 
assumptions, the analysis does not include a radial flow 
model. 


9 . 1 Expression for Qj^ in Terms of Section Loading Characteristics 


The starting points for the derivation of are equations 
(4.37) and (4.31) which are 



(9.1) 


(9.2) 



.7 
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Equations (6.63) to (6.65) define FI to F 3 in terms of 
section loads. To obtain Fi^ to from (6.63) to (6.65) 
we employ the following results for the section loads. 


V. ‘ 0 

(9.3) 


(9.4) 

ir~ 

(9.5) 

w * 

a A 

' /O 

(9.6) 

Q ^ 

(9.7) 

p A 

'Z/o ^ 

(9.8) 


(9.10) 

^ ~3 ru 

(9.11) 

* 

■ Oo) 6 

(9.12) 

(9.13) 

^0 ^ j Ut > ^ 

(9.14) 

1 T 

(9.15) 

" a ' c [4 ^ ^ 


f~ z ^ -t <j? 

(9.16) 

& ~ 0 -i- 0 ^ 

(9.17) 

B - Oc -h 0^ 

(9.18) 
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- f 


(9.19) 


z ^ > 

) 

(9.20) 

U 

- f J 

1 

a 

(9.21) 

h 

- 


(9.22) 


In the above equations, 1, d, axid m .. define aerodynamic 
lift, drag, and pitching moments. Sectors 1 emd d are 
perpendicular and parallel to the relative flow vector. 

Fig. 18, which subtends an angle O^r.to the local chord. 

The corresponding direct loads Pyio^ PzlO^ expressed 
in terms of 1 and d by means of a rotation transformation 
defined in (9.4) and (9.5). Equations (9.3), (9.7) and 
(9.8) are assumptions that Pxio^» <IylO^ q*ioA are 
negligible. These assumptions are justified by the small- 
ness of blade surface viscons tractions and components of 
normal surface forces in the xio direction, as contribu- 
tions Qj^, for largely cylindrical blades. When steady 
state ^^tion aerodynamic data eu:e used, we add to the 
other terms in the pitching moment expression, (9.6), the 
quasi-static approximation to damping in pitch, md ( p 279, 

Ref (9) ). Conditions (9.14) and (9.15) simulate forward 
and reverse flows, respectively. 

Equations (9.16) to (9.22) complete the information required 
to obtain numerical values of the loads Pyio^, Pzlo^' ' 
subsequently used to obtain an expression for QjA which is 
numerically integrated. In particular, pareuneters 
M aure Icnown when the rotor state of motion is known, and 
these in turn yield from tables the C]^, C^, and Cmg/A 
occurring in ^e loading expressions. 

The flow velocity components are referred to the X 2 . 0 ' 
defined in the next section by (9.47), and these velocity 
components are evaluated at the three-quarter chord ( (y^Q)3c/4). 
The radial ooaq>onent U^xo disregard^ in our model, ana 
does not participate in the deter^natlon of the aero- 
dynamic loads. The magnitude of Up and Up is denoted U, 
and the inclination of Up -f Ut to y^o* is denoted by^. 

Fig. 18. Expressions relating these velocity components 
to the rotor motions are derived in the next section. 
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The Xe syst«a loads required in the Fi 
(6.63) to (6.65), are obtained from the 
by means of the transfc ^nation 

to F 3 expressions 
Xio system loads 


^ X ^ A )C 
5L / C?^ A / o 

(9.23) 

with (6.33) defining 

the transformation 

(9.24) 

' 


A_j9,J. A| 

v' 


-^e 

W 


C§ 


l^dcntify fre« vector Xe with P 5 or q^. Identify 
Xxo with pxQ or q,Q, and substitute (9.3f to (9.8)^ as 
convenience suggests. We obtain fron (9.23) the X^ ayzt/m 
aerodynamic loads. These eure 


-- 

-cPc^ 


A 

to 


f - ^~r p I n- P A 

■ ■5' ' ‘ ^/C ^ I 2/0 






<p 


?r, 


f\ 

'/o 


? 


•3 


V J 


i/O 






(9.25) 

(9.26) 

(9.27) 

(9.28) 

(9.29) 

(9.30) 

(9.31) 

(9.32) 
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Tho tensions required in the Fj, to F 3 expressions, (6.63) 
to (6.65), are 


/T ^ f A . 


- 0 


(9.33) 


(9.34) 


Substitute (9.25) to (9.32), and (9.3) In the Fi to F 3 
expressions, (6.63) to (6.65). We obtain 


* 6eCv'c0^w ^ e )^A ^ 


A 

Li 

J 10 


i ®<f ' A , ) ^f^:z 1^0 j Ay 


(9.36) 


5 


t\ 

to 




/o 



^ P, 


A 


- [e/e ^ "’"[(- A 6 ^) %->'•) 


; ' 


(9.37) 
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'" 3 '' 


' ( w ' <=} X ) 




[e, ;/K- A,- A , ^ ( a , - a . ) r- 


u 


(9.38) 

with Ax and A 2 dcurivod from (5.60) and (5.63). 

Neglect the underlined terms in (9.36) to (9.38). We obtain 


r-,<^-(n{^ViUYh 


(9.39) 


) o 


Substitute (9.39) in (9.2), 2 md integrate by parts, as 
required. We obtain the aerodynamic generalized forces 



(v^- 






K? ^ ‘ ' 

4V. V,. C) 
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Pv 

u 


-- p, 


A 


'5.' 




p 4 





i 


Co 


9 


A 



^ + Fz,'t 



(9.41) 


(9.42) 


In the above, subscript T denotes a quantity at the tip, 
r » r*p. Substitution of (9.40) in (9.1) yields the aero^ 
dynamic generalized force, Qj^. In the computer program, 
the integrals in (9.40) are evaluated numerically by the 
trapezium rule. 

As a convenience, we list in (9.41) and (9.42) the loads 
required in (9.40), as well as loads csnployed to display 
the aerodynataic loads in the blade-oriented Xio system. 
The Xio system loads Pyio^ PzlO^ derive from (9.4) 
and (9.5), or from the^ transformation inverse to (9.23). 
The additional variables required in (9.42) are listed in 
(9.10) to (9.22). 

A slight var ration expresses the loads in (9.42) in terms 
of Cl, Cd, U, Up and U 7 and may be a convenience. 
Substituting in*^(9.42) the results from (9.10) to (9.21), 
as required, we obtain 

'R35 = \ 


m 
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9*2 Expressions for Relative Flow Velocity 

In this section, we derive expressions for the relative 
flow velocity at a blade section in terms of the generalized 
coordinates in physical space and their time-derivatives, 
and the velocities induced by the rotor hub motion. These 
quantities are defined when the modal displacements and 
velocities are found from the integration of the blade 
modal differential equations and the solution to the support 
or body equation. With the relative flow velocity known, 
we are aUale to define the effective angle of attack o(.^ of 
a section, (9.16), and local Mach number, M, (9.22), which 
yields the section load from tables. 

For the sake of completeness amd as a convenience, we 
include in this Section and Section 9.4, the description 
of a gust model which is preferable to that supplied 
previously under contract to NASA emd described in Ref (2) . 
This theoretical gust model was not Incorporated in the 
present program. 

We express the velocity of the flow relative to a point 
on the blade as 


U 


~v 


V 




(9.44) 


Here v is the velocity of a point on the blade induced by 
the motion of the rotor blade and hub relative to the 
stationary axis Xj. The term v&ir isolates all contribu- 
tions to u whose soxirces are independent of the rotor state 
during the calculation. Typically, v&jr may derive from a 
wind tunnel, gust, or rotor variable inflow source. In our 
scheme, the rotor inflow is Independent of the rotor state 
during the calculation and the inflow velocity is 
conveniently located in v^xr. Fig^ 19 illustrates the 
construction of tne flow velocity u for a typical case 
comprising a wind tunnel source a rotor induced 

Inflow, V* And a flow component -7 induced by the rotor 
motion. 

As a preliminary to finding U, we define the X^o Axis 
(Fig. 18) by the rotation transformation 
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Figure 19. Ex^u^>le ofjConstruction of Relative Airflow 
Velocity, U. 


217 

PAGE 


SA29 REV 0 



I 


Sikorsky Pircraft 


u 

OtVtSlON OF UNITVO AMCAAFT COMFOAATlON 



REPORT NO. SER-50912 





with 



I o o 


0 


c- 



e -6 


(9.45) 


( 9 . 46 ) 


Here € is the angle subtended by the chord and the line 
Yio' parallel to the X 5 “Y 5 plane, and lying in the Yio”*io 
plam (Fig. 18). A study of Figi^e 8 shows that € is 
approximated by the total angles, defined in (9.17). To 
preserve the rigor of our derivations, however, we do not 
assume that £ *■ 0 , and prove in the next section that our 
approximations consistently yield this result. 

Defining the component of relative airflow velocity in the 
axis as 

o'"" ^X,'» , 

with Up ^uld U>i< illustrated in Fig. 18, our objective 
reduces to finding Ut and Up in terms of our generalized 
coordinates, and the hub motions. 

Our approach will be to refer to the X 5 axis all the 
contributions to the relative air flow appeeuring on the 
right-side of (9.44), and to resolve these components to 
the X^Q* axis by the transformation 

X ,0 >^-&e (9.48) 


Substituting in the transformation A_ ^ the result of the 
next section that 


€ 



(9.49) 


and multiplying the transformation matrices in (9.48), we 
find 
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^- X ^^-6 


\ 

I 

- V 

— Vv 


V 

I 


w 

o 


(9.50) 


O 


\ 


after neglecting second and higher order products of elastic 
variables. 

Considering first the component, V^tr in (9.44), express 
this as components v^i, Vyj, V 21 referred to the Xx axis, 
and a rotor Inflow component V 21 = -V referred to the X, 
axis. The components v^j, Vyj, and V 21 are known functions 
of the spatial variables xj, yj and zj, and the time, t, 
and represent a known wind flow or assumed three-dimensional 
gust. The variable Inflow Is modeled as a known fvinctlon 
of blade radius and azimuth. 


For ^e grounded support shaft 2 ucls system of Chapter 5 , 
our Xi system velocities are 



For l^e rigid body In free flight shaft 
our Xi systoa velocities are 

(VA\'a')'2,| 


axes of Chapter 3, 


-- Tr 


-y 

(9.52) 


219 

PAGE 


1 


* “a ry n 



Sikorsky Pircraft 


u 

AMICRA^ 

R. 


REPORT NO SER-50912 



The 

transformation matrices 

in (9 

.51) and (9. 

.52) are 



f\ 


A 

- 





•r 

s 


- Ps 






^ Ps 



1 

^&s 


i 


Hs 


c 

Ps 

' ‘^'©S Vs , 





^i’s' 

\ 

c, 

As 

-f 5 (j)^ Jy/y ) 


^0$ ^fs ' 

^ Vs ^‘fs. 



(9.53) 




e J'jr 


A 


4 


‘ f\ 


-Q 


A 




Jrr 


C <j^ ‘ ^ f 


+ Ce' Sy. 
■f 


I - ' Si| 4 ^ * ' ^6 

/ 5 c^ ' 5 e ' ^ '^© 

, -C<^‘59' ^ c^'C/ 


(9.54) 

Substitute (9.53) In (9.51) and represent the airspeed in 
a wind tunnel as Vxi * Vj^lR* obtain for the grounded 
support 
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^ e ^ ^ (3.55) 

CyA»R)x , = ( ^ / s ’^^>5' ^ ^^s ^ ''^s ^ I ^ ^ 


Substitute (9.54) in (9.52) arid restrict the gust model 
to vertical gusts. We obtain for the rotor coupled to a 
rigid body in free flight 





=: -s/Ce'i^2x 

i C^'Cg'U', -V 


(9.56) 


This model allows arbitrary cylindrical gusts in steady 
horizontal translation to be represented. The gust velocity 
profile Vjx is defined as a function of distance perpendicular 
to the moving wave front. Section 9.4 describes the proce- 
dure for obtaining from this gust function the velocity Vzi 
induced by the gust at any point on an arbitrary blade of 
the multi-blade rotor. 

The program supplied under contract does not represent the 
gust in the manner described here and in Section 9.4. Details 
on the type of gust used are given in Ref. (2). 

To find the Xe components of flow velocity Induced by wind 
tunnel or gust sources, we define 
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(9.57) 

L I, Z, I 


with the right-side components derived from (9.55) for the 
grounded support mode and from (9.56) for the rigid body 
in free flightjnode, respectively. We obtain the X 5 axis 
components of from 

/ ^ij f 

(9.58) 

^ 3j (^'r^ >n')j 

with Equation (7.29) defining the direction cosines a^j in 
(9.58). This completes the derivation of Vj^xR* 

Consider ix>w the velocity of a point on the blade induced 
by the motion of the blade, we see that this follows from 
(7.2), which is 




V/ 


^ V' 


OT 


f 




Jl 


(9.59) 


Equation (7.18) yields the component^of the velocity of the 
origin of X 5 system resolved to the X^^ system. These 
components are 
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t U) 




CJ 




3/ 



'^'■'^1 *■ ^z, ^ X, ^3, ■ ‘^3 ,/^r , 


(9.60) 


with /*xl» their time derivatives defined in 

(7.24). Figure 1$ may be used to recall the physical 
meanings of the other parameters in (9.60). 


We define 

, 0 , 

C c ^ a, 3 


(9.61) 


Then by a similar derivation to that used to obtain (9.58), 
we find 




X.,' " 


^5- 

From (7.40) we have 




(9.62) 


= a 




and differentiating this we obtain 


r 


(9.63) 
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dP 


V f J-.. 


dt ‘ iOt c£t ■ ' cUt 


+ A^-r. 1% t /<■,-') /?, i/dlR,- ) P 

«t- -"i- hsT/h^' (9.64) 

Employ in (9.64), equations (7.42) to (7.44) for di 5 /dt, 
djj/dt, and dj^ 5 /dt. We obtain 

4^ r -‘"'2.,- R, , 'l 

3/^5 J 5^ 5 3^5 y 

^ J ( Pg + C02 ) 

*■ (yz .- + CO 


*T ■ ^5" 


W (9. 


65) 


with Equation (7.59) yielding ^x5' 5' andc«J25 
in (9.65). Neglect second and higher order products of 
elastic deflections in (7.63) and (7.64) and specify the 
blade section three-quarter point as the point at which 
the relative air velocity is to be evaluated, that is 




lo 



(9.66) 


Z 


( o 


- o 


We obtain from (7.63) and (7.64) 


74 0- 

r- A 

?7^’5i'>)5c/4 


V' -h 



vvP 4- 

Jc/4 

Px,- = 



PV' ^ 

V' 4 


- 

W 4 
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Tha naglact of aacond and highar order products of elastic 
variables in (9.67) and (9.68) requires the approximations 
defined by (7.72) to (7.74) to quantities appearing in 
(9.67) and (9.68). The required elements are 



■- Co. 

V/ f Se ^ w 

^ a 

‘ Ce. 

\/ + ^ vV 

% 

- 5e, 

V ^ e e 


- 

C6c 

n 

CM 

= Ce^ 

' ^e, c 

• 

- -e. 

, 5©. 

CD 

-- ^-4 

t-&e. 


Substitution of (9.59) amd <9.65)_in (9.44) yields the 
airflow velocity referred to the X 5 axis. This is 

^^5" ^ ^ Pz 5 - - ^'2 f 

U-C3. - 

(9.70) 


Employ (9.47) to (9.50) to obtain from Ux 5 » Uys, and 0^5 
the Xio' components of airflow velocity. We rind that 
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/o 




■=( 


V/ 


^rr 




- \v 


o ) X 


(9.7X) 


Up 


- w 


Substitute (9.70) in (9.71) and neglect second and higher 
order products of elastic variables. We obtain 

'L^'<)5);r,- ■ (^AIr)xs- + *lj!r 

- KVoy - (Kmr:)z^- + ^ ‘^'<5' 


p C”-' 

<^X,- --r 




(9.72) 

(9.73) 


p. 




u 


= A 








The magnitude of the flow velocity, airfoil angle of attack, 
and flow Mach number derive from 


u ■= ( U,f + uJ^) 'I? 


(9.74) 
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li - u ! 


(9.75) 

(9.76) 


Equations (9.72), and the other equations defining tt^ie 
elements in (9.72) given i\ this section and Section 9.4 
for the gust, provide the working forms for the evriuation of 
U<Ti and Up. It is to be under sl.o<^ that blade modal dis- 
placements, qj,, ard velocities, are Jcnown. These are 
employed in modal sums, like (7.80) and (7.81), to obtain 
blade physical displacements and the required spatial and 
temporal derivatives of these displacements. The s\ib- 
stitution in (9.72) of hub motion velocities ^tnd shaft 
angle displacements deriving from the grounded support or 
rigid body support system of equations completes the 
determination of Ut and up, (9.72). The parameters, U, ^j., 
and M required to calculate ^e aerodynamic loads in Qj^ 
follow from (9.74) to (9.76). 


9.3 Proof That Rotation Angle g » 6 

We prove that the rotation angle € relating and x^g 

is equal to the total angle B defined in (9.l7]r. Multiply- 
ing the rotation matric^::: in (9.48) and neglecting second 
and higher order product of elastic variables, we find 


/A 


- 6 






- Ce 

- 5, A - X 

(- 2 V 


\/ 

C (T lO- 


W (9-77) 

~ .0, 

-0 


The scalar product of the unit vectors and k^ is the 

direction cosine -b f-$ , that is 
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Since f ' is perpendicular to ke> by construction, their 
scalar product is zero, and, henSe, from (9.78) 


e - O 


(9.79) 


9.4 Derivation of Gust Inducec* Flow Velocity v^ t From the 
Prescribed Gust Function - 

We describe here the coordinate transformations for finding 
from the input gust function the velocity vzl induced by 
the gust at any point on an arbitrary blade of the multi- 
blade rotor. 

Our procedure consists of obtaining the stationary axis 
coordinates xi, yj corresponding to an arbitrary point on 
any blade of the rotor in terms of blade radial and 
azimuthal coordinates. Using these values of xj and yj 
we interpolate the gust function to obtain the velocity 
V 21 at the point r, y/ on the blade. We neglect in the 
determination of Vzi the effect on xj and yj of blade 
elastic displacements and blade pitch. The velocity tr 2 ms- 
formations of Section 9.2 then yield the contributions to 
UT and Up from Vzx» thereby completing the numerical 
definitions of U^ auid Up. 

Fig. 20 illustrates the geometry of the gust model. The 
gust is cylindrical vertical gust living with a steady 
velocity Vq perpendicular to the wave front. The gust 
vertical velocity profile in the direction zt is f . ( 
in a coordinate system, ^n, shown in Fig. 20, which moves 
with the gust, and with /xi perpendicular to the gust front. 
The length of the gust is/xi ^ instant t 3 tp 

that the gust first i^inges on the rotor disk, the origli 
0 of the moving eucis Xt^ attached to the rotor is_assumed to 
coincide with the ozigxn of the stationary axis Xx« At this 
instant, the gust is inclined at an angle or to the rotor axis 
Xx' as shown in Fig. 20, and the rotor configuration is such 
that the azimuth angle of a reference blade used as a 
time parameter, is ygt « helicopter is assumed 

to be moving always with a velocity V in the negative xx* 
direction. We wish to find xx and yx for a point 1 ^, r on 
any blade of the rotor, and thereby to obtain Vzx from f(/xl)< 
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Fig. 20 yields the transformation between the moving ^11. 
axis and the station 2 u:y Xt axis. We require only 
this is 






(9.80) 


The transformations of Section 5 yield 


(9.81? 


with XjQ_r ^presenting the Xjjcoordinates of the origin 0 
of the Xx system. Since Xjo » 0 at t » tQ we have 


^^0 

- r Vx 

(9.82) 


^ o 

- 'Vf-t -±o) 

(9.83) 


- V u 

(9.84) 



(9.85) 


= O 

Symbols v^q, VyQ ( and v^ot not required here) are 
components of nub velocity, referred to the Xj system. 


By the reversal law for transposed products, we can prove 
that 

-- Cc « -e ' in- ^ ' A e ' y I 

230 

PAGE 


run 



f 


I 




i 


Sikorsky Pircraft 


u 


REPORT NO. SER-50912 


Equations (9.86) and (9.87) enable us to use the transposes 
of (9.54) and (7.29) to obtain the left-sides of (9.86) and 
(9.87) required in (9.81). Recall from (7.28) that 
direction cosines a^^j are given by 




(9.88) 



with i in aj,j denoting a row, and j denoting a column. 
Define the new set of direction cosines bj^j from 


Lj J ■ ')fr P' 

-Ss ‘ s^'Cf t- V > 

C<p.' SQ'Cf'-t 


)rr 


- (V 5^' 

5 <p I S0> 

- C^' 5©' 


Sy.' +5^ 'Cy' 


(9.89) 

> -5^‘ 


i 


/ 


(9.90) 


with i in bj^j denoting a row and j denoting a column. We 
may express zbe components of (9.81) as 


t-j_ - -'/(•t -t„) l,,, X-, + ba, >5, 4- a , 

'iT -- \c,^ +■ baa '3, bja 




(9.92) 


Replacing t by i'^/M.axA by we summarise the 

formulas required to form 
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- Vg 
SL 

( 

■fv). ) 


-- -JT ) 4 '°n 


^ 2 / 

2; 

} / 



7 
^ 1 



^ I ^ 

r e Cy, 0 „ 




'' 











(9.93) 


Direction cosines and b^j derive from (7.29) and (9.89). 
The argument of the direction cosines aj^j and functions c^ 
and s^in (9.93) is the angle to the point r on a blade 
at which V 21 is required. It is not to be replaced by 
unless the V 21 at the reference blade is needed. 

A restriction of the gust model that should be kept in mind 
is that the rotor hub is assumed to be translating with a 
steady velocity Vxo * ~V in the calculation of the gust- 
induced velocity, V 2 j. This restriction is introduced in 
the expressions for the hub translations xio, and yjo, (9.82) 
and (9.84). The effect of the restriction is an error in 
the calculation of V 21 when the hub motion is non-uniform, 
and Vxof and VyQ exhibit time-dependence. It was thought 
that the restriction to steady hub translation in the 
calculation of the gust profile, V 21 , was worth the 
simplifications achieved in the xjo and yxo expressions, 
(9.83) and (9.85). Further, it was felt that our gust model 
is an adequate representation for problems in which the 
major translation is approximated by v^q -V. It should be 
noted that there are no hub motion restrictions in the 
velocity transformation yielding Up and u^ from V 21 , which 
include hub angular displacements for general displacements 
without restriction. Also, the use of v^o « -V to calculate 
V 21 should not be taken to mean that • -V applies to 
other parts of the theory, such as in ue inertial and other 
aerodynamic elements of the model where no such assumption 
is made. 
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A minor ^int to note is that for substantial rotor disk 
inclinations our gust will not penetrate the rotor disk 
until some time after t «■ to# and this should be understood 
to interpret a response to this gust model. 
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10. Lag Damper Generalized Force Q j 

Equations (4.38) and (4.35) are the basis for the formation 
of the generalized force induced hy concentrated loads. 

The equations are 









P. 


3 


J 


( 10 . 1 ) 



( 10 . 2 ) 


In the present model, we assume that the lag damper is the 
only contributor to Other physically concentrated 

loads, like those indiiced by the pushroa, are embodied in 
the root springs in ^e calculation of normal modes 
(Chapter 3\) , Fig. 21 illustrates the geometry assumed, 
which is such that the lag damper transmits to the blade 
only a mcnnent M 24 in the X 4 axis, that is 





r Oj O -) 



With 


(10.3) 


(10.4) 
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Figure 21. Lag Damper and Hinge Geometry. 
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Kl- - 




^ ^* 0 ' 



(10.5) 


The iTK»nent8 applied to the yoke A derive from 


'/ ( 10 . 6 ) 

from which we find 

(10.7) 


A torsion moment Mx 5 cannot transfer from yoke A to cuff B 
across the feathering bearing. It follows that Mxs » 0 at B 
and the only non- zero moment applied by the lag damper to 
the blade at B is 




Z 





( 10 . 8 ) 


The moment per unit length of blade follows from ( 6 . 68 ), 
which yields 






(10.9) 


To obtain this result we assumed the existence of the 
derivatives of M 25 on the grounds that the applied loads 
eure distributed physically over a small but non-vanishing 
region of blade, in such a way as to make the derivatives 
of M 25 finite. Substituting (10.8) in (10.9) and assuming 
M 25 ' to be predominantly larger than Mz5, we find 


t 



( 10 . 10 ) 
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Substitution of (6.63) to (6.65) for 
yields 




J ' ^ 




REPORT NO. SER-50912 
to ?3 in (10.2) 


( 10 . 11 ) 


( 10 . 12 ) 


(10.13) 


Assume that all derivatives in (10.13) are finite, 
integrate (10.13) by parts, and set v^j = 0 at r = 0 
and qzs' » 0 at r ■ tt* Wo find 



where ^ is a small distance over which is applied the moment 
per unit span, qz 5 » induced by the lag deunper. We approxi- 
mate this moment as 

r C y 

/3 " 

Approximating the blade deflections as 


(10.15) 

(10.16) 
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w 




(10.17) 


(10.18) 


vdiere subscript zero Indicates a quantity at the lag hinge, 
substituting (10.17) and (10.18) in (10.11) and (10.14) and 
taking the Units as A-^ 0, we find 


C- 



(10.19) 

( 10 . 20 ) 

( 10 . 21 ) 


Assume the lag damping moment representation to be 

( 10 . 22 ) 

(10.23) 



where Cxj) is the lag damper constant (lb~ ft- sec/rad units) , 
and substitute (10.23) in (10.21). We obtain 


0 




'o 9a 


r- 

O 


(10.24) 
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Substitution of (10.19), (10.20), and 
yields 


which completes the definition of QjP 
element tj, defined by (4.77). 
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11. Shears and Moments 


We derive here the expressions used in the program to obtain 
the running moments along the blade and the force and moments 
applied by the rotor to the hub. The running moments are 
employed to display the blade response. The rotor forces 
and moments applied to the hub are employed to calculate 
the excitation of the support to find the response of the 
coupled rotor/support system. The coupled system includes 
either the grounded flexible support or the coupled fuse- 
lage in free flight. 


Blade running moments for display of program responses are 
calculated from the following approximate internal moment 
expressions. 



r 

(11.1) 



(11.2) 


_/o ^ 

'A ' 

(11.3) 



(11.4) 

- <iJl^ R, R, 


(11.5) 


Expressions Ill.l) to (11.3) derive from (6.17), (6.30) and 
(6.31) with T siibstltutlng for F^io* certain conditions 

- for example, when hub accelerations are important - t 
should be replaced by Fxig (See (6.61) and (6.71) to obtain 
FjjIo) • Also, the expression for q i« not strictly correct 
when a Si)cors)cy-type counterweight is present, and preferably 
should be replaced by (6.31). 

Blade root shears and moments resolved to the axis derive 
from the external force and moment expressions, (6.71) and 
(6.74) to (6.76), specialized for x > 0. The resulting 
shears and moments are subscripted with the symbol e. 
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c ^ 


( 11 . 6 ) 




(11.7) 


^ )c(i^ V JJ(v-^(6:))(p/^.u%-'')c(i< 

- S^Yw-w{'»))(p>55 + 
Ja/-'xx(.'>) /(P x? *■ fjfj ) 

-S\x-o)(fA^P-^s)^ 


The Inertia loads Px5^ to qzsl^ in (11.6) and (11.7) are 
obtained from (8.14) to (8.17) with (7.7?) yielding ^ 
Cz/ and (5.63) and (5.60) yielding \ andA 2 . Aerodyneuaic 
loads Px5^ to qz 5 ^ are obtained from (9.25) to (9.32) wiUi 
(9.4) to (9.6) yielding Pyio^» Pzl0^» *nd qxlO^ required 
in these expressions. For articulated rotors the program 
sets the hinge moments to zero and loops around the calcu- 
lation of (11.7). 
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The corresponding shears and moments resolved to the axis 
are 


h' 

). 





(11. 

8} 


)? . 



■V 




(fzi 

). 






and 

this 

may be 

written 







('pA,-)e 

■ '^311 1 


P “"i/ (Pzy )e 

(11. 

9) 

CP'S,)? - 




^ ^9 Je 




- 


( )e. 

P ‘^Aj( 


p '*'}} 



with 

i the 

Sj^j obtained from (7.29). 

Similaurly 



(P'/c,)e 


CMx,-) 

e ^ 

(",sV 

F %j(ht^-\e 

(11, 

,10) 



,K,-) 






1 a. 

\ V 


)e T 

.4- yC 




Fig. 22 may be used to derive hu' shears and mome.. s induced 
by a single blade. The rcsult.3 i 
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( 11 . 11 ) 


(M„)h = (M.,). + ^ 'f 


Shears and moments from all the blades eure summed to obtain 
forces and moments applied to the hub by the rotor. 

The forces and moments in (11.11) and (11.12)_ are components 
in the direction of the shaft oriented axis Xx used with the 
grounded support (Fig. 4) . To obtain components in the 
directions of the shaft oriented axis f i coupled to the 
rigid body (Fig. S) required for the determination of body 
response in free flight, we apply the following additional 
transformations . 



■ >F^, ,-F.)^ 


(11.13) 


(11.14) 


These relations may be verified from (5.12) or Fig. 6. 
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An alternative to (11.7) is to express the moments at the 
hinge as internal reactions consisting of the moments 
illustrated in Fig. (23) . Resolve these to the X 5 eocis. 
We obtain for the moments applied by a blade to the rotcc 



(11.15) 

which reduces to 

I - ^Lo S' ^ (U.16) 

C^2^)e - S 


The motivation for using the external forces to calcula 
the shears from ( 11 . 6 ) was a desire to increase the acc ccy 
of the root shear calculation in comparison with the internal 
force method. The latter method would require third deriva- 
tives of mode shapes and for accurate results usually many 
more blade modes than are adequate for the response calcula- 
tion would be needed. On the other hand, the uce of (11.16) 
instead of (11.7) for the hinge moments would be more 
efficient computationally and should be accurate, except 
possibly for very stiff springs. Also, the more accurate 
external running moment expressions, (6.74) to (6.76), 
should replace the Internal running moment expression, (11.1) 
to (11.3), now used in the progreun. 
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Figure 2 . Internal Reaction Moments at the i»iade Root. 
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12. ExpresBion for Tim»-D«pand«nt Conponyit of Blade Pitch 
Angle/ and Non-Diaeniilonalization Rules 

We derive here the expression used for the time- dependent 
component of blade pitch angle, and the rules employed 
in the program to non-dimensional ize the modal equations. 

12.1 Pitch Angle Component, & y 

The expression for the time-dependent contribution to the 
pitch angle induced by cyclic control inputs emd pitch- flap 
and pitch-lag couplings is 

' 9 3 (12.1) 

Angles Ais and are cyclic inputs. The remaining terms 
^are derived below, and the version of (12.1) used in the 
program is defined. 

Fig. (24) illustrates the geometry assumed to obtain the 
pitch- flap coupling. Point A is the attachment of the pitch 
horn. Point P is the attachment of the pushrod. If the 
point P is unrestrained and free to move with the blade, 
it moves to P^., following a flapping deflection and elastic 
bending deflection at A. A rigid pushrod restrains P such 
that it stays at P 2 , and this geometric configuration is 
assumed . 

The change in pitch angle is 


A9 

A 




(12.2) 







^ (24 

)pl 

- P4! 


(12.3) 


) ^ 




(12.4) 


'P, - 




(12.5) 
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Substitute (12.3) to (12.5) in (12.2). We obtain 

Af - ''p)) (12.6) 


p 


Expand the elastic deflection as a Taylor series from the 
blade hinge. We obtain 








(12.7) 


( 12 . 8 ) 


From (3.12) and (3.14) we find 










(12.9) 

Approximate ^by (12.9) amd substitute (12.7) to 

(12.9) in (12.6). Add to these expressions for the 
cyclic input contribution and the program input pitch lag 
coefficient tanci^^^. We obtain 


( 12 . 10 ) 


^ t ~ fs^ ^r' ' ^ S’ , - S 


^ ^ V. oC 

- s 


3 
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'f 
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■h 


pi':>p 

S' - 


ei, 


Ct-I 


( 12 . 11 ) 




O^-v^ o( 


r ^ 

o( -)- — 








o 


where tanoCi^ is the input pitch-lag coefficient. 

In the program the pitch horn is assumed to be located at 
the hinge (^X = 0) . The elastic deflection, is assumed 
to apply at the mid-point of the first segment of the blade. 
For these assumptions, (12.10) reduces to 

~ ' ^ / 5 ^ 'f - ■ /^ 5'j S ~tiK^ zU 

- ^ f\ 5 ^ 

(12.12) 


12.2 Non-Dimensionalization Rules 

The purpose of this section is to define the rules for 
obtaining non-dimensional forms of the variables appearing 
in the blade modal equation. These variables eure used in 
the coding of the modal equations. With these non- 
dimensional variables, the modal equation in the non- 
dimensional space have exactly the same forms as the 
dimensional modal equations, cited in the earlier text and 
no corrections reflecting changes in form are required. 

To obtain a non-dimensional variable, we select a combina- 
tion of terms from among iOq,!!, and R which will non- 
dlmensionalize the variable. Here mo is a reference mass 
per unit length of the blade (slug/ft) . We replace 
derivatives with respect to time, d/dt, by d/dy^. We 
replace radial derivative d/dr by d/d(r/R) . Frequ(*.ncies 
are replaced by frequency ratios, a>/Hand rotor speed /I is 
replaced byH/Aa i. we denote variables non-dimenslonalized 
in this way by overbars, and use an arrow to denote a sub- 
stitution. The following are examples of this transformation. 
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Appendix 14.1 - Coefficients in Expression for Correction to Modal Stiffness, 
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Appendix Ik. 2 - Derivatives of Direction Cosines a 
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